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ABSTRACT 

We present an overview of the HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) m the 
Magellanic Clouds project* \\ Inch is a Herschel Space Observatory? open lime key program. We mapped the Large 
Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) at 100, 160, 250. 350. and 500 /an with lire Spectral 
and Photometric Imaging Receiver (SPIRE) and Photodetector Array Camera and Spectrometer (FACS) instruments 
on board Herschel using the SPIRE/PACS parallel mode. The overriding science goal of HERITAGE is to study 
the life cycle of matter as traced by dust in the EMC and SMC. The far-infrared and submilUraeter emission is an 
effective tracer of the interstellar medium (ISM) dust, the most deeply embedded young stellar objects (YSOs), and 
the dust ejected by the most massive stars. We describe in detail the data processing, particularly for the PACS data, 
which required some custom steps because of the large angular extent of a single observational unit and overall 
the large amount of data to be processed as an ensemble. We report total global fluxes for the LMC and SMC and 
demonstrate their agreement with measurements by prior missions. The H LRU AGE maps of the LMC and SMC 
are dominated by the ISM dust emission and bear most resemblance to die tracers of ISM gas rather than the stellar 
content of the galaxies. We describe the point source extraction processing and the criteria used to establish a catalog 
for each waveband for the HERITAGE program. The 250 //m band is the most sensitive and the source catalogs 
for diis baud have ^25*000 objects for the LMC and ^5500 objects for die SMC. These data enable studies of 
ISM dust properties* submillimeter excess dust emission, dust-to-gas ratio. Class 0 YSO candidates, dusty massive 
evolved stars* supernova remnants (including SN1987A), Hit regions, and dust evolution in the LMC and SMC. 
All images and catalogs arc delivered to the Herschel Science Center as part of the community support aspects of 
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die project. These HERITAGE images and catalogs provide an excellent basis for future research and follow up 
widi other facilities. 

Key words: catalogs - dust, extinction - infrared: galaxies - Magellanic Clouds - submillimeter; general - surveys 
Online-only material: color figures 


1. INTRODUCTION 

The Large Magellanic Cloud (LMC) and Small Magellanic 
Cloud (SMC) are excellent aslrophvsical laboratories for study- 
ing the life cycle of the interstellar medium (ISM). Located at 
50 and 60 kpc, their proximity permits detailed studies of re- 
solved ISM clouds and their relation to stellar populations on 
global scales. In tire Milky Way (MW), which is difficult to study 
in its entirety, ISM studies are more challenging because they 
suffer Ironi line-of- sight crowding and distance ambiguity that 
requires extensive work to disentangle (e.g., Reid et al. 2009). 
In the LMC and SMC, ISM clouds and stars are at a similar 
distance, rendering luminosities directly comparable. In partic- 
ular, the thin-disk morphology and favorable viewing angle of 
the LMC make the comparison between dust emission and gas 
line emission fairly straightforward (Table 1 ). The Magellanic 
Bridge is a filament of neutral hydrogen, which joins the SMC 
and LMC over some 1 5 kpc (Muller et al. 2003a; 1 Imdman et al. 
1 963; Mathewson 1 985; Staveley -Smith el al. 1 998). lire part 
of the Magellanic Bridge that is closest to the SMC has been 
called the Tail, a term that we adopt for this paper. The Tail has 
evidence of recent star formation (Harris 2007; Grondin et al. 

1 992; Hambly et al. 1 994; Chen et al. 20 1 3). When we refer to 
the mapping observations of the SMC, we are referring to die 
SMC+Tail region 

The LMC and SMC differ substantially in many global 
properties such as in eta Hi city, molecular gas content, and 
gas-to-dust ratio thus allowing us to study how die processes 
governing galaxy evolution depend on these properties ( fable I ). 
T he physical processes leading to star-formation (in particular, 
dynamics, H-H 2 transition, gas cooling) can be examined at 
metallicities of 0.5 (LMC) and 0.2 (SMC; Table 1) which 
is similar to galaxies at the qioch of peak star formation 
in the universe (Madau et al. 1996; Pei et al. 1999). The 
lower nieta 11 icily results in lower dust content in both galaxies 
compared to the MW (Table 1 ). Less dust in the ISM increases 
the mean free path for stellar photons from the massive 013 star 
population which in turn increases the ambient UV radiation 
lieldsto higher values than in the solar neighborhood. In terms of 
galaxy evolution, the LMC -SMC pair is well suited to studying 
how the agents of evolution, die ISM and stars, operate as a 
whole in two galaxies dial are tidal ly interacting with each odier 
and w idi die MW (e.g., Zaritsky & Harris 20 04; Bekki & Chiba 
2005; Besla ct al. 2007; Gardiner & Noguchi 1996). 

Full surveys of the LMC and SMC provide both a detailed 
view of aslrophvsical processes occurring within these galaxies 
and an opportunity for statistical assessment of these processes 
over an entire galaxy. Moreover, a whole-galaxy survey pro- 
vides a complete template for high-redsliift galaxies. The SMC 
and LMC have been surveyed al many wavelengths revealing 
structures on all scales thereby providing a rich context for inves- 
tigations. The stellar content has been photometrically surveyed 
at optical (MCPS; Zaritsky el al. 2002, 2004) and near- infrared 
(near-IR) wavelengths (Two Micron All Sky Survey (2 MASS), 
Skrutskie et al. 2006, Nikolaev & Weinberg 2000; Vista Magel- 


lanic Cloud Survey, Cioni el al. 201 \ ) and for stellar variability 
(MACHO, Alcock et al. 2000; IRSF, Ita et al. 2004; OGLE III, 
Udalski et al. 2008a, 2008b). The ISM gas content has 
been traced in atomic hydrogen widi the Hi 21 cm line 
(ATCA+Parkes, Rim ct al. 2003; Parkes, Slanimirovic ct al. 

[ 99 Muller ct al. 2003a), in ionized gas content in die SI IASS A 
survey by Gaustad et al. (2001), and in molecular gas content 
in the CO rotational transitions in die NANTEN, MAGMA, and 
mini -south Columbia surveys (Fukiri el al. 2008; Wong el al. 
20 1 1 ; Mizuuo et al. 200 1 ; Cohen et al. 1 988: Rubio et al. 1 99 1 ) 
and spectroscopically for 1 1 1 and 1 1 2 surveys (Tumlinson el al. 
2002). Molecular and atomic gas fractions vary greatly between 
the LMC and SMC ( fable 3 ). The dust content, relative to the 
gas. lias been analyzed in UV spectra along lines of sight to 
bright stars (Gordon ct al. 2003; Maiz Apellaniz & Rubio 20 12). 
The emission from dust has been surveyed with IRAS (Schwer- 
iug Sl Israel 1 989), and the Midcourse Space Experiment (Egan 
et al. 200 1 ). The extinction by dust has been studied by Dobashi 
et al. (2008), 

In our Spitzer Surveying the Agents of Galaxy Evolution 
(SAGE) surveys, we have invested >1000 lir fully surveying the 
LMC and SMC over the full wavelength range (3.5-160 ^m; 
Meixner et al. 2006; Gordon et al. 2011) with die Spitzer 
[RAC and MIPS cameras and with die Infrared Spectrograph for 
hundreds of point sources and for tens of taigeted ISM regions 
to study the lifecycle of galaxies (Kemper et af 2010). The 
Spitzer Survey of the SMC (Bolatto et al. 2007) pre- dated the 
SAGL surveys and covered the brightest regions of the SMC. 
These Spitzer studies were very successful in providing galaxy- 
wide overviews of the emission of warm dust in the ISM (e.g., 
Bernard et al. 2008). discovery of thousands of young stellar 
objects ( YSOs; e.g. : Whitney et al. 2008; Gruetidl & Chu 2009; 
Carlson et al. 2012; Sewilo et al. 2013), and a census of the 
mass injected by lens of thousands of dusty evolved stars (e.g., 
Srinivasan et al. 2009; Boyer et al. 2011, 2012; Riebel et al. 
2012). Hie analysis of die dust emission in die ISM compared 
to die gas emission reveals higher dust-lo-gas ratios than for 
the UV extinction measurements (fable 1 ) suggesting eidier 
missing gas mass as traced by H 1 2 1 cm or CO or there are large 
variations in the dust-to-gas ratio. 

!n this paper, w e present an overview of the HERschel Inven- 
tory of The Agents of Galaxy Evol ution ( HERITAGE ) program, 
an open time key program on the Herschel Space Observa- 
tory * ( Pi 1 brad et al. 2010 ). HERITAGE provides an imaging- 
photometric survey with the Photodetector Array Camera and 
Spectrometer (FACS; Poglitsch el al. 2010) and Spectral and 
Photometric Imaging Receiver (SPIRE; Gridin et al. 2010) at 
100, 1 60, 250, 350, and 500 pm. These crucial long wavelengths 
complement die Spitzer Space Telescope studies ol' die LMC 
and SMC. We discuss the observing strategy in Section 2, while 
die data processing steps are described in Sections 3 and 4, 
Sections 5 and 6 describe the spatial distribution oft he imaged 
emission and the Magellanic Clouds' total integrated Herschel 
emission, respectively. T he source catalog creation is described 
in Section 1, and die makeup of the catalog is discussed in 
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Table I 

Probities qf the Large Magellanic Cloud and Small Magellanic Cloud 


Physical Properly 

LMC 

SMC + Tail 

References 

Distance (kpe) 

50 

61 

1,2 

Disk thickness (kpe) 

2.5 

4.2 

3 

Inclination 

23' -3 7* 

62* 

3 

Metall icily (Z Q ) 

0.5 

01-0,2 

4,5 

ISM gas size (kpe) 

8 

3 

6, 7. 8 

Hi gas mass (Mg) 

4.8 x 10 s 

4.2 x 10 s 

9, 7, a 

CO gas mass (A/g) 

4 x 10 7 

4 x 10 s 

10,11,12 

Dust-to-gas ratio (lA) 

0.005 

0.001 

13 

Dust-to-gas ratio (IR/radio) 

0.0017 

0.0001 1-0.00083 

14. 15. 16 

Stellar mass (Mg) 

L7 x tO 9 

3.7 x 10 s 

1 7, 18 

Star formation rale (Mg yr l ) 

0,1-0,25 

0.024-0.05 

17, 1 9. 20, 21,22 


References. (1) Ngeow & Kan bur 2008; (2) Szcwezyk el al. 2000; (3) Submit ani an & Subramaniam 2009, 2010, 
2012; (4) Russell & Dopila 1992, corrected by Asplund el al. 2004; (5) Lee el al. 2005; (6) Kim el al, 3 998; 
(7) Slanimirovie el al. 1999; (8) Muller et al. 2003a; (9) Staveley-Smilh el al. 2003: (10) Eukd el al. 2008; (11 ) 
Mizuno et al. 2001; (12) Muller el al. 2003b; (13) Gordon el al. 2003; (14) Bernard el al, 2008: (15) Bol et al, 
2004; (16) Gordon el al. 2009; (17) Skibba et al. 2012; (18) Harris & Zarilsky 2004, 2009; (19) Whitney el al, 
2008; (20) Sewilo el al. 2013; (21) Bolatto el al. 2011; (22) Wilke cl al. 2004. 


Table 2 

Principal Characteristics of die I HERITAGE Survey, Herschel Program ID KPQT mmeixner 1 


Parameter 

LMC 

SMC 

R.A. (J2000) 

Q5 h 18 n, 48 li 

01 h 30 m l6 s 8 

Deck (J2000) 

6834T2" 

-73 J 22'58'.'8 

Areal coverage 

8° x 8? 5 

5° x 5° + 4° x 3° 

ADR sizes 

30' 36' x 450' 495' 

80' 102' x 120' 240' 

Epoeh 1 dates 

2010 Apr 28 lo May 2 

2010 Mar 24 25 

Epoch 1 OBSID^ 

1342195668, 1342195669, 1342195683 

1342192680. 1 34219268 L 1342192697, 


1342195684. 1342195707, 1342195708 

1342192698. 1342192699 


1342195712, 1342195713, 1342195728 


Epoch 2 dales 

2010 Aug 2-7 

2010 Jim 18-21 

E|>odi 2 OB SID 

1342202086, 1342202087. 1342202202 

1342198565, 1342198566. 1342198590, 


1342202203, 1342202216, 1342202217 

1342198591, 1342198863 


1342202224, 1342202225, 1342202243 


Epoch 3 b dates 
Epoch 3 Of) SID 

1342202244 

2010 Sep 23-24 
1342205049, 1342205050, 
1342205054, 1342205055, 1342205092 


Notes. 

a OBS I D is the reference number for (tie AGR executed and enn be used lo identify I he da la in die Herschel science 
archive (HSA). For epoch 1. we nole that for the SMC, these OBSID only include (he PACS dala (see text). 

'■ Wc have an epoch 3 dala set for the SMC only because it was a makeup for I he epoch 1 . Epoch 3 has bolh PACS and 
SPIRE dala. 


Section 8. We conclude with a description of the higher level 
dala products available to the community in Section 9 and a 
summary in Section 1 0. 

2. OBSERVATIONAL PLAN 

Using 285 hr of Hers chei observations, we have performed 
a uniform photometric survey of the LMC (8 x 8? 5), SMC 
(5 g x 5° and 4 g x 3 g ) in all SP1RL bands (250, 350, 500 / 2 m) 
and in tire 1 00 and 1 60 /mi PACS bands. The raw data observa- 
tions have been stored in the Herschel archive under proposal 
KPOT_mmeixner_l where KPOT stands for i£ key program in 
open time” The map sizes for die TMC and SMC were cho- 
sen lo fully encompass each galaxy as observed at 160 /mi 
(SpitzerfhALPS, Figure 1 ). Observing characteristics for die sur- 
vey have been outlined in fable 2. 

For lire full SP1RL/PACS imaging oflhe LMC and SMC, our 
science goals were primarily driven by the need to spatially 


map the coldest (T < 15 K) dust which requires SPIRE, 
The ability to acquire high-quality PACS images at 100 and 
1 60 /mi for a minimum of extra time using the SPIRR/PACS 
parallel mode fuddled secondary science drivers getting better 
wavelength coverage ( 100 /mi) and spatial resolution ( 1 60 /mi ) 
oflhe large, warm (25 K < 7’ < 1 5 K) grains and compact source 
spectral energy distributions (SLDs) than prior observations 
by IRAS and S/V/m/MIPS. We efficiently and fully mapped 
the LMC and SMC at 100, 160, 250, 350, and 500 /tm using 
the SP1 R F /PAC S para Helm ode at low sc atm i ng speed ( 2 0 " s 1 ). 
lire low scanning speed mode was picked as it allowed us to 
reach our sensitivity goals, to distinguish between glitches and 
sources in die timeline, and to produce minimally smeared PACS 
point spread functions (PSFs). 

the astronomical observing requests (AORs), ranging from a 
few to 9 hr in duration, were designed to scan over each galaxy in 
long strips to ensure that background measurements were made 
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Figure I, HERITAGE SPIRT AOR mapping strategy tor the LMC (top) and 
SMC (bottom) overlaid on the MIPS 160/^m image from the SAGE-LMC 
Spitser l egacy program (Meixnei el a! t 2006) and SAGE- SMC Spitzer Legacy 
program (Gordon et a!. 201 1 ). Each long red rectangle represents the footprint 
of one AOR (see Table 2 for details on sizes). The long scans of each epoch 
arc sufficient to gel off (lie main galaxy far-infrared emissions and sample a 
suitable background. Hie AORs for bolh epoch 1 and epoch 2 arc shown. 
Because the epochs were separated in lime by approximately 3 months, the 
AORs are orthogonal and the combination of both epoch tool prints results in a 
basket weave appearance. This observing strategy improved the linat photometry 
quality of the combined-epoch maps. 

(A color version of litis ligure is available in the online journal.) 


at the beginning and end of each scan leg, allowing removal 
of the 1 // noise with timescales longer than the time to take 
a single scan leg. To further suppress the 1 j f noise, each map 
was repeated with the scan angle rotated 90 ' , which naturally 
happened in about 3 months at die positions of the LMC and 
SMC. To eiTcct this orthogonal rotation in die mapping we 
scheduled die observations to have a d tree -month separation 
between epochs I and 2. This observing strategy is similar to 
the strategy that we adopted for die MIPS observations on die 
SAGE-LMC (Meixner et al. 2006) and SAGE-SMC (Gordon 
et al. 2011) projects. However, we note that it differs from 
other large-area mapping programs with Herschel. e.g., 111GAL 
(Mol ina ri el al. 2010), which used smaller AOR tiles and die 
1TSPOT “orthogonal*’ option to obtain a perpendicular direction 
as part of die AOR. Our observ ing strategy had some unforeseen 


consequences for the PACS data processing within die Herschel 
Interactive Processing Environment (HI PE) which we discuss 
below. 

The scan angle rotates by 1 day 1 for die Magellanic Clouds, 
so we concatenated several AORs in order to make a complete 
map without gaps in the minimum time al each epoch. For die 
LMC, the epoch ! date range was 20 1 0 April 28 to May 2 and 
die epoch 2 dale range was 2010 August 2-1 (Figure 1). For 
die SMC, die epoch I date range was 20 1 0 March 24 and 25; 
however* aldiough the PACS data were successfully obtained, 
die SPIRE observations failed due to commanding issues. Thus 
for the SMC, we rescheduled observations widi epoch 2 date 
range of 2010 June 18-21 mid epoch 3 date range of 2010 
September 23 and 24. Aldiough the LMC mapped region is 
approximately square in shape, we used non-uniform sized strips 
to cover the LMC in an optimal fashion (f igure 1 and fable 2) 
in order to meet our time constraints. The rectangular SMC 
coverage was broken into two approximately square shaped 
regions for the mapping die SMC and Tail regions (Figure t and 
Table 2), 


3. PACS DATA PROCESSING 
3. L Background 

Hie HERITAGE data were obtained in scan mapping mode. 
Scan maps are performed by slewing die spacecraft at a 
constant speed along parallel lines (see chapter 5.1 of die PACS 
observer's manual). In the following, we refer to a scan leg 
as one of these scan lines, i.e., a continuous segment of data 
acquired between turnarounds. During each AOR, die arrays of 
PACS bolometers (32 x 16 bolometers for the !60/nn band 
and 64 x 32 for die 100/im band) sample the far-infrared 
(FIR) dust emission as a function of time as the spacecraft 
performs die scans. As a result, raw PACS data consist of a 
three-dimensional cube, die first two dimensions corresponding 
lo die spatial dimensions of the array, and the third dimension 
to lime, which can be directly transformed into position on the 
sky using the pointing information. We refer lo a frame as the 
data in the PACS bolometer array at lime \ (a frame is therefore 
two-dimensional, 2D), and to a timeline as die data recorded 
in one bolometer as a function of time (a timeline is one- 
dimensional, ID). At the end of the data processing, the 
timelines therefore need to be projected onto the sky, leading 

10 the creation of a 2L> map of the MR emission recorded by 
PACS. 

The PACS bolometers are affected by the so-called 1 // 
noise caused by thermal drifts where die amplitude of the 
noise varied inversely with frequency. Measurements of the 
FACS low frequency noise indicated it approximately followed 
a I j yfj relation (Sauvage & Okumura 2009). In addition, the 

11 LRU AGE AORs were the longest in angular extent for the 
Herschel mission, with an extent of about ^8' in length for the 
LMC, and -5° for the SMC. Asa result, the HERITAGE PACS 
timelines are dominated by baseline drifts resulting from 1// 
noise, which therefore need to be mitigated prior to pro jection of 
die data onto the sky. In addition to diese drifts, cosmic ray hits 
on bolometers and readout electronics cause glitches or sudden 
jumps in the signal strength of the timeline. These glitches also 
need to be removed during the data processing. 

Noise, drifts, and glitches can be estimated and removed using 
the redundancy of die observations, since each region of the sky 
is seen by several bolometers during the scans and cross- semis. 
Software programs based on the redundancy are available to 
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mitigate drifts and l // noise, such as MADmap (Cantalupo 
et al. 2010), Scanamorphos (Roussel 2012) and Tamasis (P. 
Chanial et al, in preparation). Because they are based on 
redundancy, these algoridims require PACS data to be processed 
in its entirety in one run. In order to use such algorithms on 
the large amount of data collected by the HERITAGE survey 
would have required substantial computing resources, perhaps 
a supercomputer. For example the LMC PACS 1 00 /xm data 
takes 200 Gbyte, which was well beyond our 1 28 Gbyte RAM 
server For a future data release, we will work with the Tamasis 
program on a supercomputer to generate PACS images of (lie 
LMC and SMC. However, for the initial data release we have 
processed the data using our own custom pipeline that included 
some steps in addition to the standard PACS data processing 
software in order to appropriately remove the baseline. Our 
custom approach fits within the computational resources of our 
128 Gbyte RAM server Our initial approach creates images 
which liave residual l // noise that limits the scientific utility 
of the data for some projects, e g., power-spectral analysis of 
the ISM structure. TTie different steps of this custom algorithm 
arc described in the following subsections. Figure 2 illustrates 
the steps of the PACS data processing, especially the baseline 
subtraction, for one timeline corresponding to one bolometer. 

3.2. Basic Reduction /torn Level 0 to Level 1 

The basic PACS data reduction was performed starting from 
the level 0 product to the level 1 product using the standard 
pipeline software in die HIPE version 7.0 data reduction soft- 
ware (Ott 2010; PACS data reduction guide - fi ). First, calibration 
blocks are identified in the telemetry and processed. Rad and 
saturated pixels are masked. Readouts arc converted into volt- 
ages. Pointing information and sky coordinates are added to die 
frames. The signal (in volts) is converted to Jy using die respon- 
sivity factors in 1 11 PE 7.0, and a flat field correction is applied. 
After these steps, we exported the data as FITS files for further 
processing within the IDL software environment because die 
HIPE pipeline software was unable to accomplish the task. 

3.3. First-order Baseline Subtraction 

The level 1 frames in die timeline data were compromised 
by 1// noise which was caused by drifts in the bolometer 
signal w ith time. T here were also large pixel -to-pLxel offsets 
between bolometers. This made the removal of glitches, which 
were sudden jumps in the timelines due to cosmic rays hilling 
the detectors, very inefficient widi the algorithms implemented 
in HIPE. As a result, a first-order, rough baseline subtraction 
was implemented (in IDT) before applying the deglitching 
algorithm. For this purpose, we split the signal timeline of each 
bolometer by scan leg. For each scan leg and each bolometer, 
we estimated the median signal at the end points of die scan, 
located outside of die LMC or SMC, where we assumed that 
the FIR Hux was zero. Hie median of die signal at die end 
points was estimated in a 600 frame long ( 1 0 f ) window tor 
the PACS 160 jum band, and 300 frame long (10 ) window 
for the PACS 1 00 /mi band. Regions of the timelines located 
w ithin 1 000 frames of a calibration block w ere masked out to 
eliminate the influence of transients, which were drifts in die 
signal following a calibration block. A linear lit to die median 
values of the timelines at die scan cuds was then estimated, and 
subtracted from the timeline corresponding to the scan leg of 


The PACS mev guide is on line at hi Ip: herschd. esac.esa.int 4icsSrdoc-7.0/. 


each bolometer. This process is illustrated in the top two panels 
of Figure 2, w here the red line indicates the linear fit to the end 
points of each scan leg. This rough baseline subtraction also 
removed pixel- to -pixel offsets. 

3.4. First Step of Deglitching 

Alter this first-order baseline subtraction was applied, the 
multiresolution median transform (MMf) deglitching algorithm 
was implemented to remove the effects of cosmic ray hits, Tlie 
details of the MMT algorithm implemented in 111 PE have been 
described in die PACS data reduction guide and references 
dierein. We used die recommended parameters for the PACS 
1 60 /an band, with mcates = 2, nsigma = 5, incr-faci = 2, 
mode = multiply , For the PACS 100 /mi band, we observed 
that applying those parameters led to the core of the PSF of 
bright point sources also being removed. I fence, we selected a 
more conservative set of parameters for the green band, with 
n scales = 1 . This approach preserved more of die source flux 
at die expense of more glitches remaining in die PACS 100 pm 
images. The residual glitches in the 1 00 pm band dierefore need 
to be removed at a later stage using the second level deglitching 
algorithm. 

3.5. Refined Baseline Subtraction 

The next step in the PACS data processing is to remove die 
contribution of iff noise and baseline drifts more accurately 
than in the rough baseline subtraction. We developed our own 
method in IDL to remove iff noise and drifts, which takes 
advantage of previous, lower resolution FIR observations of the 
Magellanic Clouds with IRAS at 1 00 pm (Schwering 1989), 
Spitzer MIPS at 1 60 /mi (Meixner el al. 2006; Gordon et al. 
20] l) and DIRBE at 140 pm and 240 pm (Silverberg et al 
1993 ). The basic idea of die method is that the baseline of 
the PACS timelines can be estimated as the difference between 
the PACS timelines smoothed to the resolution of the MIPS 
or IRAS observations and a synthetic MIPS or IRAS timeline 
created from the coordinate timeline of die PACS frames and 
the MIPS or IRAS images. Here and below when we refer to 
die MIPS 160 /im images, we mean a composite MIPS 160 pm 
and DIRBE image. The areal coverage of die PACS 1 60 pm 
image is larger than the MIPS 160 pm image. We extended 
the MIPS 160 pm image widi an interpolated llux distribution 
derived from the DIRBE 140 pm and 240 pm images. Note dial 
small differences in the central wavelengths between the PACS 
1 00 and IR IS 100, or PACS 160 and MIPS 1 60, combined with 
color corrections are negligible (of order 2% of typical dust 
temperatures) and dierefore cause only small discrepancies in 
die baseline subtraction. 

In order to prepare synthetic PACS- 1 ike timelines, the MIPS 
160 pm and IRAS 100 pm images of the LMC or SMC arc first 
converted from MJy sr _1 to Jy pixel 1 , where die size of the 
pixel corresponds to the size of the PACS 160 /an bolometer 
(6:4), or of the PACS 1 00 pm bolometer (372). A MIPS or IMS 
synthetic time hue is dien created as the value of die converted 
MIPS or IRAS image at the R.A./decl. location in each PACS 
160 /an frame of die timeline. The PACS 160 pm and PACS 
1 00 pm timelines are then convolved using a 1 D convolution to 
MIPS 160 pm and IRAS 1 00 pm resolutions, respectively. l r or 
the PACS 160 pm band convolution, we use a slice through the 
center of the convolution kernel created by Aniano et al. (20 1 i ), 
while for the PACS 100 pm band convolution, we use a slice 
through the center of a Gaussian convolution kernel of FWHM 
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Figure 2. FACS timelines for one bolometer and during one AOR, at the different steps of the FACS processing algorithm. Stalling from the lop. the first panel shows 
the raw timeline in black, on which is over-plotted in red the linear fit to the end points of each scan leg. performed during the rough baseline estimate. The second 
panel shows the timeline resulting from this rough baseline subtraction. l"he third panel shows the timeline alter initial rough baseline subtraction, smoothed to MIPS 
resolution using a ID slice through the center of the 2D convolution kernel in black, ihc synthetic MIPS timeline obtained from the FACS mapping scheme and the 
MIPS map is show n in red. The fourth panel shows the MIPS timeline in red. and in black, She FACS timeline, after a rough baseline subtraction, smoothed to MIPS 
resolution, and after performing a linear interpolation inside masked bright and/or compact regions. The lifih panel shows the final estimation of the residual baseline 
(he,, after rough baseline removal from a linear fit to the end |>oints of each scan leg), obtained by subtracting the MIPS timeline to the smoothed, interpolated FACS 
timeline in the fourth panel. The final, baseline subtracted FACS timeline is shown in the bottom panel. 

(A color version of this figure is available in the online journal.) 


width v^FWHM 7 ff> , s — FWHMf ACS ), where FWHM /jms =4.'3 
is the FWITM resolution of the IRAS 1 00 /an image* and 
I WllMpAcs — 8" is the resolution of the PACS 100 /an image. 
An example of smoothed PACS timelines is shown in the third 
panel of Figure 2. 

Smoothing the timelines in 1 13 was not equivalent to extract- 
ing a 1 D timeline from a 2D smoothed map but was necessary 
to boot-strap our way to the hrst 2D PACS maps. Nevertheless, 
our approach caused some artifacts in regions with bright and/or 


compact sources that required some special attention. Because 
the convolution occurred in I D, the signal “dilution” is not as 
high as in the 2D convolution case. As a result, the shape of 
bright and/or compact sources in the smoothed PACS timelines 
did not match the shape of the same sources in the Ml PS 1 60 f t m 
or III IS 100 /an limelines. This effect was manifest in Figure 3, 
which showed a zoom on a point source in PACS 160 /an and 
MIPS 160 /an timelines. Because the baseline was estimated 
as die difference between the smoothed PACS timelines and 
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Figure 3. Baseline subtraction and point sources. Top: example of point source 
morphology mismatch between the smoothed PACS and the MIPS 160 jum 
timelines. The black line shows ihe roughly baseline subtracted PACS timeline, 
corresponding to the timeline in the second panel in Figure 2* The red line 
shows the MIPS timeline, as in the third panel of Figure 2. The blue line 
shows Ihe roughly baseline subtracted, smoothed PACS timeline, as panel 3 
of Figure 2. Bottom: effects in ihe projected maps of over-subtraction around 
bright /compact sources due to differences in ihe shapes ol the MIPS and PACS 
timelines caused by the one-dimensionality of the convolution. This problem is 
solved by li nearly interpolating the timeline inside such regions. 

(A color version of this figure is available in the online journal.) 


ihe MIPS 160 /mi or IRAS 1 00 pm timelines, this mismatch 
therefore resulted in over- subtraction lobes on either side of a 
compact or bright source, as shown in Figure 3. Maps of the 
subtracted baseline showed that this problem also affected the 
flux of more extended sources in a subtler way. To avoid this 
issue, we masked point, compact, and bright extended regions 
in the baseline subtraction process. The masks are shown in 
Figure 6, as an example for die LMC. To estimate the baseline 
inside those regions, the smoodied FACS and the MIPS 160 pm 
or IRAS 100 fim timelines were linearly mteqiolated between 


die signal in die frames located right outside a given region, and 
the baseline was estimated as tire difference between the inter- 
polated, smoothed PACS and MIPS 160 /mi or IRAS 100 pm 
timelines. Outside the masked regions, the baseline was simply 
estimated as the difference between the smoothed FACS time- 
lines and die MIPS 160 /mi or IRAS 100 pm timelines. This 
process was illustrated in the fourth panel of Figure 2. The 
baseline was subtracted from the PACS timelines. An example 
of baseline subtracted PACS timeline is shown in the bottom 
panel of Figure 2. 

3.6. Removal of Long Glitches in the PACS / 00 pm Band 

We found the presence of discontinuities in the brightness 
images for the 100 /an band (see Figure 3). These discontinuities 
were caused by jumps in the timelines affected by cosmic rays 
and affect an entire line of the PACS 1 00 pm detector. These 
so-called long glitches were not removed by die deg li telling 
algorithms implemented in HIPE (MMT and second level 
deglitching). We therefore developed a custom IDT algorithm 
to remove these artifacts. First, the rough position of the glitches 
was estimated by recording their coordinates in die PACS 
images. Second, the total signal in die timeline of the detector 
was computed, and each timeline around the coordinates of 
die jump was sc aimed by eye. The shape of die jump in die 
total timeline was obvious, as shown in Figure 6, allowing a 
detenu i nation of the affected line's location (a.k.a. index). To 
detenu ine the index of die frame at w hich the glitch occurred, we 
computed the difference between the timelines and the timelines 
shifted by ! and 2 frames. The mdex ol die jump was computed 
as die index where dris difference is greater than 0.2 Jy pixel 1 . 
The sections oldie tun dines on the I ell and right side of the long 
glitch, which show a characteristic ramp caused by the cosmic 
ray hit, were lit to a polynomial. T he fitting was done iteratively 
in order to try different orders of polynomial and lengths around 
the glitch to perform the lilting. Once a reasonable order and 
length for the tit w ere identified, typically a poly nomial of order 
two and a lit 150 frames long on either side of the glitch, the lit 
was removed from the timeline of each bolometer in die affected 
detector line around the glitch. The result ol tins algorithm was 
shown in die right panel of Figure 3. 

3. 7. Correction of A stnmiethc Offsets 

At the time of the observations, the Herschel spacecraft still 
suffered from degraded pointing reconstruction performance. 
This degraded pointing caused astrometric offsets of up to 
5"-l 0" in the PACS images, which is significantly larger than 
the nominal astrometric accuracy of Herschel is 2 !f at the la 
level 3 These astrometric offsets posed problems not only in 
the combination ol die different epochs available (e.g., smearing 
of die PSP, double sources, etc.), but also in the combination 
of the HERITAGE images with previous Spirzer images of the 
Magellanic Clouds (Meixner et al. 2006; Gordon et al. 201 1). 
Therefore, w e developed an algorithm (in ID!,) to correct those 
astrometric offsets based on die positions of point sources in 
die MIPS 24 pm catalog, the astrometric accuracy of which is 
tied to 2 MASS and accurate to MX' 3 (Sknitskie et al 2006). We 
chose to use die 160 pm band of PACS radier Ilian its 100 pm 
counterpart because the 1 60 pm band is deeper and because the 
baseline removal is superior for that band, being based on the 
MIPS 160 map at higher resolution than the IRAS 100 map. 


liilp: hcrschd.esac.esa.inl Docs FACS html ch03s06.htm] 
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Figure 4. Spiizer MIPS 1 60 /j. m (lop) ji n cl IRA S 100 /cm (bottom) maps of the 
LMC. The masks Used to interpolate Ihe PACS timelines inside bright or compact 
sources in Ihe LMC during the baseline estimate are overlaid as blue circles. 
The 50 MJy sr -1 surface brightness contour, corresponding to Ihe transition lo 
the nonlinear regime of MIPS 160 /tm is shown in green in the MIPS image. 


First, cadi scan leg of each AOR in (lie PACS 160 /nn data 
was mapped to spatial coordinates in R.A. and decl. using the 
PhotProject algorithm in I UPF. PSF fitting was performed on 
each scan leg image using the model PSF, and a PACS 160 /an 
point source list was subsequently extracted. For each scan leg, 
the PACS 160 /xm point source list was cross-correlated with 
the list ol point sources identified in the MIPS 24 /xm catalog, 
llie (constant) astrometric olTset for each scan leg was then 
calculated as the median difference between the point source 
coordinates in the PACS 160 /an images and in the MIPS 24 gin 
images. The astrometric olfset was finally subtracted from the 
coordinates of the frames. Asa final check on the astrometric 



Figure 5. Example of a long gliteh. a cosmic ray liil aCTccling an enlire line of 
bolometers in Hie PACS 100 /xm delcclor, ihe lop panel shows Ihe image of 
the glitch, while ihe bottom panel shows Ihe projected limelines after the long 
glitch has been removed by our custom algori thm. 


1.0 



- 1.0 

1.018* 10 1 1.020x1 0 s 1 .022x10® 1.024x10 s 1.026x10 s 

Frame Number 

Figure 6. Manifestation of a long glitch in the PACS 100 ^m timelines. The 
black line shows the uneorrected timeline, with the discontinuity at frame 
number 1 .021 x 10 s caused by the glitch. The red line corresponds to a 

polynomial Eil of order two to Ihe 150 frames on each side of Ihe discontinuity. 
Ihe blue line shows Ihe corrected timeline (glileh removed). 


correction, w e extracted sources from the PACS 1 00 and 1 60 /n n 
images (Section 7 below') and did another comparison lo Ihe 
MIPS 24 //m source locations. We found a small 
but systematic offset between the PACS sources and the MIPS 
24 fim sources hi the R.A. and decl. directions that we corrected 
for both the PACS 100 and 160 /mi images. Figure 7 shows the 
histogram of the source offsets between PACS 160 and MIPS 
24 found in the LMC and SMC after die correction. The final 
astrometry of the PACS images is consistent both with the MIPS 
24 u m and by extension the SAGE surveys, and with the SPIRL 
data sets. 


3.8. Mapping ami Second Level Deglitching 

As pointed out in Section 3,4 S residual glitches subsisted in 
the PACS 100 /on images, Iliose were eliminated by applying 
the second level deglitching algorithm in 111PE. We used the 
sigma-clipping algorithm in median mode with a box size of 
10 and a sigma threshold of 3 The PACS 160 /xm images were 
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Figure 7. Left panel: the histogram ofposiiional oITsels between the PACS 160/im point source locations and their corresponding MIPS 24 /jm point source locations 
for the LMC. Right panel: similar positional histogram offsets but for the SMC. 


Table 3 

Principal Characteristics oltlie HERITAGE Survey; Herschel Program ID KPOT mmeixner 1 : Wavelengths, Resolution, and Sensitivity 


Characteristic 

PACE A (/tm) 


SPIRE k (/i ni) 


100 

160 

250 

350 

500 

Instrumental pixel size (") 

3.2 

6.4 

6 

10 

14 

Image pixel size (") 

2 

3 

6 

10 

14 

Angular resolution, FWHM ("), RA C) 

6.7 x 6.9, 62.3 

10,7 x 12.1,93 

!E,2 

24.9 

36.3 

Sensitivity,, 1 (kr (MJy sr 1 ) 

85-90 3 

32-50“ 

6 

3 

2 


Note. a llte noise in PACE varies in the map with masked regions having the higher noise values Ilian the unmasked regions, see text 
(Section 3. 10) for details. 


successfully deg Inched using the MMT algorithm, so this slop 
was unnecessary for the 160 /nn band. The deglitched, baseline 
subtracted, and astrometrically corrected PACS 160 /an and 
1 00 fim timelines were mapped using the PhotProject task in 
HIPF. Maps for the LMC were created from the combination 
of both epoch 1 and 2 data. For the SMC, maps were created 
by combining the data from epochs 1-3. The limi ting surface 
brightnesses listed in Table 3 represent 10tr of die standard 
deviations measured on the outer regions of the images where 
there is no significant signal. For both the LMC and SMC, we 
used the individual epoch maps to under stand die uncertainties 
for the source extractions from the PACS images. 

3.9. Cmss-calibralion of PACS 100 and 
160 tim with Prior Missions 

We examined the accuracy of the absolute flux calibration 
with PACS. For this purpose, we convolved and resampled die 
PACS 160 /nil (or PACS 100 /mi) images to the same resolution 
and grid as the MIPS 1 60 /an (or IRAS 1 00 /an) images. For die 
convolution to a common resolution, we used the convolution 
kernels provided by Aniano et al. (20 1 1 ) for the PACS 160 /an 
images, and a Gaussian kernel of FWHM 43 (IRAS 100 /an 
resolution) lor the PACS 100 /an images. We evaluated and 
applied color corrections for PACS. MIPS, and I HAS, and found 
respectively 0.97, 0.99, 0.97, and 0.96 for PACS 100, FACS 1 60, 
MIPS 160, and IRAS 100 respectively, assuming a modi lied 
blackbody with temperature 25 K and emissivity spectral 


index = 2. We then performed a pixel-to-pixel correlation 
between the PACS and MIPS /IRAS llux to cheek for gain 
differences. 

We found that PACS and MIPS 1 60 /an agree within 4% 
for surface brightnesses below 50 MJy sr 1 . The top panel of 
Figure 8 shows the pixel-to-pixel correlation between the PACS 
1 60 /an and MIPS 160 /an maps, A bisector linear lit to the 
correlation using only the data below 50 MJy sr -1 yielded a 
slope of 0,96 (i.e., PACS = 0,96 MIPS), shown with a green 
line in Figure 8 . Note, that, without applying any color cor- 
rection, this s!oj>e becomes 1 .00. The agreement is expected 
since ( 1 ) the baseline subtraction algorithm effectively sets the 
calibration of the PACS maps to the MIPS calibration out- 
side of the regions used to mask bright and/or compact FIR 
emission in the PACS timelines (see Section 3.5), and (2) 
those regions roughly correspond to the 50 MJy sr" 1 contour 
of the MIPS 160 map (Figure 4). Note that color correc- 
tions were not applied during tire baseline subtraction pro- 
cess, which explains why the slope of the MIPS 160 /PACS 
160 pixel-to-pixel correlation is closer to ! when no color cor- 
rections are applied. Although the match between the PACS 
160 and MIPS 160 flux calibration below 50 MJy sr 1 is 
built-in our baseline subtraction technique, Paladini et al 
however showed in a report on the PACS extended flux c a li- 
bra! ion (lit tps: / / nhsc sc i . i pac , cal tech . edu/pacs/ doc s/ Pho lometer/ 
PICC-NI lSC-TR-034.pdf) based on PACS maps processed with 
Scanatnorphos or MADtnap. that PACS 1 60 and MIPS 1 60 were 
in agreement within 5%-20% below 50 MJy sr -1 . 
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Figure 8. Pi xe l-l 0 -pixel surface brightness correlation between the combined epochs of the PACS 160 //in convolved, resampled images of the l.MC, and! the MIPS 
160 /jtn on the one hand ( top), and the interpolation at 160 ft m between the D1RRE 140 and 240 /urn bands on the other hand (bottom). Surface brightnesses are plotted 
in MJy sr 1 . In [he top panel, the green line shows a bisector linear lit, of slope 0.96. lo the correlation in the linear range of the MIPS detector (surface brightness 
below 50 MJy sr -1 ). flic red line shows an exponenlial fit lo (lie eonctalioTi in the nonlinear range of MIPS (surface brightness above 50 MJy sr =l ), which yields 
(PACS/50 MJy sr“ 1 ) == 0.96 (MIPS/50 MJy sr -1 ) , J 1 in the nonlinear regime. In Ihe bottom panel, (he red line shows a I ; 1 correlation. 

(A color vei'sion of this figure is available in the online journal.) 


Above a brightness threshold of 50 MJy sr 1 s which corre- 
sponds to regions inside the masks described in Section 3,5 S 
there was a departure from a linear correlation between MIPS 
and PACS 160 /fin; (lie MIPS 160 /fin brightness was lower 
than tire PACS 160 /fm brightness. Recall dial, because die 
baseline was linearly interpolated within die masks marking 
the 50 MJy sr -1 contour of MIPS 1 60 and was not based on 
the MIPS synthetic timelines, a match between the MIPS 160 
and PACS 160 is not necessarily expected in this case. The 
flux calibration inside those regions is left to the PACS calibra- 
tion, The PACS response function is known to be linear up to 
1G 4 MJy sr" 1 , while MIPS is linear only up lo 50 MJy sr -1 . Thus 
this departure between die PACS and MIPS flux calibrations 


above 50 MJy sr s was most likely due to die nonlinearity of the 
MIPS detector. A power- law fit between the PACS and MIPS 
160 /an brightnesses for the flux ranges above 50 MJy sr 1 
gave (PACS/50 MJy sr" 1 ) = 0.96 (MIPS/50 MJy sr" 1 ) 1 n (red 
line in Figure H). Because die flux calibration of PACS inside 
the 50 MJy sr" 1 contour of MIPS 160 is not tied to the flux 
calibration of MIPS 160, the MIPS nonlinearity was therefore 
not propagated to die PACS maps via die baseline removal 
algorithm. 

Previous COBE DIRBF observations of the LMC and SMC 
provide an additional independent way to validate die flux 
calibration of our PACS /.an maps. Although DJRBF does not 
have a 160 /an band, we estimated the brightness at 1 60 /an 
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Figure 9. Pixel-lo-pixel correlation between Lire convolved, resampled PACS 100//m combined map of the LMC and the IRAS 100/vm image on the one hand (top) 
and I he DIRJ1F, 100 //m map on the other band (bottom). Surface brightnesses are plotted in MJy sr '.In the top panel, the green line shows a bisector linear fit to the 
correlation, yielding a slope of 0.97. The red line shows an exponential fit, leading to (PACS/50 MJy sr 1 ) ~ 0.94 (IRAS/5 0 MJy sr 1 ) 1,os < In the bottom panel, ihe 
red line indicates a 1 :1 con elation. 

(A color version of this figure is available in the online journal.) 


by performing a linear interpolation in Ihe log-brightness of 
the 140 and 240 /cm bands. After applying the proper color 
corrections, we analyzed the pixel-to-pixel correlation between 
the D1RBF interpolation at 160 /mi and the PACS 160 /mi 
brightness, which is shown in the bottom panel of Figure 8. 
The red line shows a 1;! correlation. Despite the large scatter, 
presumably caused by the peculiar sampling scheme of the 
DIR 131: maps, the flux calibration of PACS 160 /mi appeared 
to be consistent with the flux calibration of DIR RE across the 
whole brightness range. This is an independent confirmation 
that onr PACS maps are well calibrated, since the baseline 
removal was in no way tied to DIR RE. In addition, COBH 
was calibrated on mi extended source (the cosmic microwave 


background (CMB)) and therefore provides the best extended 
calibration source. 

We performed a similar analysis on the PACS 100 /an 
band. The pixel-to-pixel correlation between the IRAS 1 00 /an 
and PACS 100 /an maps is shown in the lop panel of 
Figure 9. There was good agreement between the PACS and 
IRAS calibration, A bisector linear lit to die correlation over the 
whole flux range yielded PACS = 0.97 IRAS 2.18 MJy sr -1 . 
An exponential fit over die whole flux range yielded 
(PACS/50 MJy sr” 1 ) = 0.94 (IRAS/ 50 MJy sr 1 ) 1 *® So 
there was a slight nonlinearity in this pixel-to-pixel correlation. 
We also examined die pixel-to-pixel correlation between the 
DIRRE 1 00 /an maps, and the correspondingly re-convolved. 
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Figure 10. Zoom-in on the PACS 160 {lop) and 1 00 (bottom) jum maps of the 
LMC. The blue circles correspond to masks used to estimate the 1// noise and 
bolometer drills (Section 3.5, Figure 4). The green box shows an example of 
region used to estimate the residual noise in masked regions of the LMC, 

resampled PACS maps (see bottom panel of Figure 9), The 
pixel-to-pixel correlation between the liual PACS 100 ptm maps 
and (lie DIRBE 100 fim was consistent with a 1:1 correlation. 

3,10. Noise Estimation 

Because the 1 / / noise is estimated and removed in different 
ways inside and outside the masks define by the 50 MJy sr 1 
contour of MIPS 1 60, residual noise in the PACS maps is not 
homogeneous. We have estimated the residual noise in those 
different eases by computing the standard deviation of the 
PACS brightness in several regions, inside masks but outside of 
obviously visible ISM structure on tire one hand (see example 
shown as a green box in Figure 10)* and outside masks in the 
outskirts of the Magellanic Clouds on the other hand. For both 
the LMC and SMC, we found a 1 a nns of 20 MJy sr 1 
and 9 MJy sr ] in masked regions, and 3.2 MJy sr 1 and 
8,5 MJy sr 1 outside of the masks for the 160 and 100 /mi 
respectively. Residual l / / noise appears as striping in the PACS 
maps. The difference in the quality of the Iff noise removal 
between masked and unmasked regions is obvious in Figure 10: 
the striping is of larger amplitude in the masked compared to the 


unmasked regions. Outside the masks, the 1 // noise is estimated 
as the difference between PACS 1 60 (resp. PACS 1 00) timelines 
smoothed to MIPS 160 (resp. IRAS 100) resolution. Thus* the 
I // noise is removed down to size scales comparable to the 
MiPS 160 (respectively IR/iS 100) resolution for the 160 and 
1 00 pm bands. Inside the masks, the l // noise is calculated 
as the difference between linear interpolations of the smoothed 
PACS and MIPS timelines, where the points used to perform the 
linear interpolation a re located tight outside the mask. Inside the 
masked regions* the 1 / / noise is therefore corrected on length 
scales greater or equal to the radius of the mask, w hich is larger 
than die resolutions of MTPS 1 60 and IRAS 1 00 urn bands. 
Since die amplitude of 1 jf noise increases with size scales, 
the residual 1 // noise, and hence striping, inside the masked 
regions is significantly larger than outside the masked regions. 

3.1 L Limitations and Benefits 

Due to the peculiar nature of the removal of drifts and 1 / f 
noise in the 1 1ER1TAGE data, there are limitations to the possible 
applications of this product. In particular, die noise properties are 
inhomogeneous. This precludes studies of the power-spec tmm 
(or similar) of the ISM. In addition* residual striping makes the 
study of the structure of compact objects* such as 11 u regions and 
supemovae remnants, difficult, Nonedieless* we have derived a 
technique based on tittering of the Fourier spectrum to remove 
this residual striping on small compact regions. In the 2D 
Fourier spectrum of a small portion (a few aremin in side) 
of the PACS maps, the striping very clearly appears as a line 
perpendicular to the direction of die striping in die image. After 
being identified in die Fourier spectrum, the striping can be 
inverse Fourier transformed, and subtracted from the image. 
This technique has been successfully implemented in Matsuura 
et ah (201 1 ) and Otsuka et al (2010)* and does not affect the 
flux calibration. Additionally, we arc currently developing an 
algorithm, Tamasis (P. Chanial et al, in preparation), similar to 
Scanamorphos or MADmap, to remove 1 // noise in the PACS 
timelines. Processing the HERITAGE data through Tamasis 
requires significant time on a super-computer, which makes die 
testing and development of diis software a lengthy process. At 
the moment, the accuracy of die flux calibration and different 
arti tacts present m die resulting maps are being worked on. We 
will provide Tamasis maps lo the community in the future. 

Asa check on our approach, we have applied both our 1 // 
noise removal algorithm and Scanamorphos to publicly avail* 
able maps of M3 1 (taken as part of the Herschel Exploitation of 
Local Galaxy Andromeda: Fritz et al. 2012 ), and odier galax- 
ies from the KINGFISH (Key Insights on Nearby Galaxies: a 
Far-Infrared Survey with Herschel: Kennicutt et al. 201 1 ) and 
SINGS (Spitzer Infrared Nearby Galaxies Survey; Kemiicutt 
et al. 2003 ) samples, to check whether our custom algoridim 
introduced flux calibration issues via the lie to previous ob- 
servations with different detector technologies. We found the 
PACS flux calibration resulting from botli algorithms were cor- 
related widi a 1 : 1 relation* and that the PACS flux calibration of 
maps reduced with Scanamorphos, independently from Spitzer, 
IRAS * or DIRBE, was in excellent agreement ( 10%) with those 
previous missions. 

Because die baseline estimation and removal are based on 
previous FIR observations of the Magellanic Clouds w idi Spitzer 
and IR IS, potential calibration issues in Spitzer and IRAS may 
be propagated in the PACS maps. I iowever, we have verified dial 
the flux calibration of the HERIT AGE PACS maps is consistent 
widi previous measurements (COBE DIRBE, Spitzer, IRAS) 
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in regions where die baseline is estimated independently 1'rotn 
those previous observations (Section 3.5). Independent studies 
(e.g, 7 Paladini et al. 2012) have also shown dial the PACS 
flux calibration is consistent within 5%-20% with in the more 
“diffuse** brightness regime (he., not masked in our baseline 
estimation). We emphasize that die Herschel PACS observations 
constitute tremendous progress compared to the previous Spitzer 
ones, thanks to the improved spatial resolution (a factor of 3.5 
improvement for die 160 /mi baud). Our baseline subtraction 
algorithm does not affect the native PACS resolution in any way, 

4. SPIRE PROCESSING 

Hie SPIRE images of the Magellanic Clouds were reduced 
using version 7 of 11IPE, with the addition of customized 
routines to subtract the background, adjust the astrometry, 
identify glitches, remove residuals from temperature drifts, and 
mask discrepant data. 

4. L Processing Jwm Level 0.5 to Level I 

We developed an algoriOim to identify “concurrent glitches” 
he., glitches that affect all bolometers of a given SPIRE array 
simultaneously, which were most probably due to impact of 
cosmic rays on die bolometers' silicon support. We donated this 
algorithm to the SPIRE Instrument Control Center (ICC), and 
it was included in the standard pipeline. 

We observed that during strong temperature variations in 
die SPIRE cryo-cooler, die relation between die drift of die 
bolometers* signals and die plate temperature had a slightly 
different slope compared to ihe standard relation. For Uiis reason, 
the temperature drill correction implemented in die standard 
pipeline left residual drifts in die signal (Figure II), lliese 
residuals were significant in our data, especially in die 500 pm 
band, given the lengdi of the scans (^9° in the I ,MC ). We derived 
new temperature drill coefficients for scans performed during 
strong cooler temperature variations measuring the thermistor 
voltages and die bolometer signals at each end of die scans. This 
approach reduced die residual drifts considerably. 

4,2. Astrometry Correction 

We corrected the aslrometiy of SPIRE data using an approach 
similar to the one we uti lized for die PACS data ( see Section 3.7) 
except on an AOR basis as opposed to a scan leg basis. We 
used the MIPS 24 /mi source catalogs of die LMC (Meixner 
et al. 2006) and the SMC (Gordon el al. 20 1 1 ) as the common 
astrometric standard for both PACS and SPIRE because we were 
processing these data sets in parallel and needed a common 
astrometric reference frame. Moreover, the 250 /mi image 
which is our most sensitive band in die HERITAGE survey 
has many sources not found in the PACS bands because PACS 
was significantly less sensitive Ilian SPIRE. The SPIRE 250 /mi 
catalogs have many more sources in common with die 24 /mi 
catalogs, which are also very sensitive, mid thus the 24 /mi 
catalogs provides very good astrometric reference frames for 
HERITAGE. 

For each AOR, we generated a map at 250 /mi and measured 
the positions of point sources using the IDL code Stailinder 
(Diolaili et al. 2000). These positions were dien compared lo 
positions in die 24 /an catalogs of die LMC (Meixner et al. 2006) 
and the SMC (Gordon et al. 201 i ) lo detenu ine the astrometry 7 
adjustments to the data. Hie first epoch of the LMC was only 
offset by a fraction of an arcseeond from the 24 /an catalog 
in all AORs, so no adjustment was made to diose data. Hie 


science demon stration phase (SDP) observations and second 
epoch of the 1 ,MC T and bodi epochs of the SMC, displayed 
offsets of several arcseconds that differed from AOR to AOR, so 
the astrometry for those AORs was adjusted to bring the average 
positions in line with the 24 /an positions before die final map 
was produced. Ihe final distribution of offsets was not well- 
described by a Gaussian, but was peaked at Y!5 in both galaxies, 
with a long tail to offsets of several arcseconds (Figure 12). We 
measured the offsets at the longer SPIRE wavelengths to be the 
same as those al 250 /an lo widiin the uncertainties implied by 
die larger beams: so identical corrections were applied al those 
wavelengths, Asa final check on die astrometric correction, 
we extracted sources from the SPIRE 250, 350, and 500 /mi 
the combined-epoch images (Section 7 below) and did anodier 
comparison lo the MIPS 24 /an source locations. We found no 
systematic offsets which confirms that astrometry is consistent 
bodi with die MIPS 24 /an and by extension die SAGE surveys, 
and with the PACS data sets. 

4,3, Background Subtraction, De striping, and Masking 

t he background subtraction was performed separately for 
each scan hue. Hie level subtracted was defined by a linear 
function that tits a single data point on either side of the galaxy, 
each of which was the result of a sigma -clipped median after 
masking out die region of the galaxy. The mask was die same 
for all three wavelengths and was detenu hied using a simple 
map at 500 /mi. the wavelength a l which the apparent extent of 
eacli galaxy was largest. 

Since the observed regions were dominated by the extended 
emission of Magellanic ISM, we applied to die bolometer 
timelines the extended emission relative gain factors. These 
gains, detenu ined by the SPIRE ICG, represent the ratio between 
die response of each bolometer to the extended emission and 
the average response, i he difference in the response of each 
bolometer was mainly due to the variation of die beam area 
among bolometers. The use of those relative gains decreased 
the noise in die final maps by several tens of percent at high 
(>few times 1 0 MJy sr~ ] ) surface brightness. 

We then applied a destnpmg routine to background- 
subtracted data. This routine was a modified version of the 
destriping task implemented in H1PE. The routine measured the 
diff erence between die signal measured by each bolometer for 
each scan and the signal on the reconstructed map a l the same 
sky coordinates. This difference, as a function of tune along the 
scan for each bolometer was fitted with a second-order poly- 
nomial, and subtracted from the bolometer's signal timeline. A 
new map was reconstructed and the procedure was repeated for 
1 00 iterations. 

Our destriping tontine also included die detection of jumps 
in bolometer timelines and a second level deglitching. Jumps 
in bolometer timelines w ere similar to PACS long glitches (see 
Section 3.6). but in SPIRE they affected one bolometer at a time 
(since each bolometer has its own readout electronics), lliese 
jumps were detected by performing a wavelet transform with 
a Haar wavelet of fixed width of die difference between the 
signal measured by a bolometer in a scan and die signal on the 
reconstructed map at die same sky coordinates. If the maximum 
of the result was larger than a given dueshold ( spec i lied as 
quanta of the noise), die bolometer timeline was affected by a 
jump and flagged to be omitted in the map making. The number 
of flagged timelines was <0.5%. Second level deg li telling 
was implemented to remove glitches left over by the first 
level wavelet deglitcher, which worked on bolometer timelines. 
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Figure 11. Top panel: a section of the SPIRT 500 jum band map of the LMC jwoduced with the standaid 1 RPR 7 pipeline. The vertical and hori zonial stripes are due to 
residuals left hy temperature drift correction during strong temperature variations (see Section 4). Bottom panel: same region as in the top panel of the SPIRE 500 
band map produced using new tem|>erature drift coefficients. 

(A color version of 1 his figure is available in the online journal.) 


Our second level deglitching algorithm identified glitches by 
comparing die signal of a sample lo the corresponding signal 
of the reconstructed map. A sample with a glitch differed from 
the reconstructed map by more than 5 times die map pixel's 
error value comprised of the instrumental noise and the standard 
deviation of die combined samples. This algorithm did not Hag 
samples taken near or on bright sources. The number of flagged 
samples was around 1 0-20 per scan leg per array. The de striping 


routine reduced all residual drifts to less than the noise level, 
except for a single leg in die LMC, which we removed from 
subsequent processing. 

4,4. Mapping 

Map reconstruction was done with the SPIRE default map- 
making algorithm (simple map making). We converted the maps 
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Figure 12. Top panel: the histogram of positional offsets between the SPIRE 
250 /an poinl source locations and Iheir corresponding MIPS 24 fun point 
source locations for the LMC. Frequency is the number of sources with a 
particular offset. Bottom panel: similar positional histogram offsets but for 
the SMC. 

(A color version of this figure is available in the online journal.) 

output from HTPE 7 which are in units of Jy beam 1 , to MJy sr 1 , 
by dividing by beam areas (using die revised SPIRE instrument 
beams posted in 2012 December) of 1 1.093, 1.932,4.156] x 
10“ 8 sr or 1465, 822, 1768] arcsec 2 al [250, 350, 500] /an. A 
rough estimate of the uncertainly on lhe whole-galaxy measure- 
ments was determined by comparing the total llux density of 
each galaxy to interpolations between the prior submillimeter 
data of the LMC and SMC compiled by Aguirre et af (2003), 
file SPIRE data we present here were found to be consistent 
within the errors for both die LMC and SMC, respectively. In 
addition lo die two SPIRE epochs, we incorporated die 1 8 In of 
SDP data, which covered 2 J x 8 strip through the LMC center, 
at a position angle of ^22?5, described by Meixnerel al. (2010) 
into die dual LMC SPIRE map but not lor the final FACS map 
because of die bad striping of the FACS SDP data. The limiting 
surface brightnesses listed in Table 3 represent 1 Oo of die me- 
dian value of the noise maps produced during the SPIRE data 
processing. 


5. SPATIAL DISTRIBUTION OF LMC 
AND SMC HERSCHEL EMISSION 

The images resulting from the HERITAGE project provide 
the highest spatial resolution view of the emission at Ihese 


wavelengths to date. Figure 13 shows the HERITAGE data 
in all five bands for die LMC with the MIPS 70 /nn image 
from SAGF-LMC (Meixner et al 2006). At first glance, all the 
images look identical which indicates diat the emission from 
all these bands arises from die same sources of ISM dust, lhe 
images show a disk of dust emission with complex li lain ciliary 
structure punctuated by brig lit compact regions. The brightest 
of ihese compact regions, located on the eastern side, slightly 
so utli of center are the highly active star formation regions of 
30 Doradus and N160 slightly further south. The other bright 
regions are also associated witii active star formation w here the 
bright young OB stars efficiently heal die ISM dust. On closer 
inspection, one notices that die contrast between these bright 
regions and the fainter diffuse ISM emission decreases with 
increasing wavelength. This contrast becomes more apparent in 
die three-color image of the PACS 100 and 1 60 /Lin data and 
die SPIRE 250 /nn data (Figure 14). If all tlie images had an 
identical morphology on all scales, the images would appear as 
a gray scale The color differences apparent in the image ret! eel 
the differences in temperature of the dust emission. The bright 
white spots are the hottest regions, the blue or green regions are 
warm, and the red regions are coolest. 

t he l ILR1TAGL data for the SMC reveal a more disturbed 
structure w ith a NE-SW bar to the west and an extension to the 
east southeast that is called die "\ving 5? of the SMC (Figure 1 5)< 
The similarity of the PACS and SPIRE band emissions underlies 
die common source ofISM dust emission. West beyond die wing 
is the beginning of the Magellanic Band between the LMC and 
SMC, tlie so-called Tail of die SMC, and very liltle emission 
is apparent in tins region. A low r level of wispy emission is 
apparent throughout die SMC field of view and the majority 
of tliis wispy emission is attributable lo the MW's IR cirrus 
emission (e.g. 7 Gordon et al. 2009). Part of die reason for the 
low er emission in the “Tail” is undoubtedly due to a decrease in 
dust-to-gas mass ratio in die tail compared to the wing and bar 
of die SMC (Gordon et al. 2009). lhe SMC ISM dust emission 
also has a filamentary appearance with bright knots of emission, 
especially hi die bar, such as N66 in die northeast pari of the 
bar. file diree-color image (Figure 16) highlights lhe bright 
knots in white as the hottest dust regions. The Herschel spatial 
resolution is sufficient to reveal die mini-spiral shape of N66 in 
die northeast part of the bar Tlie three-color image of the SMC 
also reveals a range of hot, warm, and cool dust regions within 
die galaxy (Figure 1 6). 

Figures ! 7 and 1 8 compare the spatial distribution of the 
Herschel emission hi the LMC and SMC lo other galaxy - wide 
tracers of the stars and gas. The Spitzer I RAC 3,6 /urn band 
traces the evolved stellar population, which have SED peaks 
in die near-lR (Meixner et al 2006; Gordon et al 201 1). The 
SHASSA Her image traces the ionized atomic hydrogen gas 
(Gaustad et al. 2001 ). Hie ATCA/Farkes 1 1 1 2 1 cm image traces 
the neutral atomic hydrogen gas in the LMC (Kim et al. 2003) 
and SMC (Stanimirovic et al. 1 999). TheNANTEN CO J= 1-0 
emission traces molecular gas in die LMC (Fukui et al 2008) 
and SMC (Mizuno et al. 2001 ). 

Hie distribution of lhe Herschel emission does not resemble 
dial of the evolved stellar population, the source of the stellar 
wind and dust production, which is traced by' the I RAC 3.6 /mi 
band. Indeed the LMC bar, w hich is dominated by light from the 
evolved stars, is not apparent in the Herschel images. Herschel 
emission has a similar morphology to ISM gas tracers, especially 
1 1 1 . indicating that Herschel emission arises primarily from 
the cool dust of die ISM. However, die dynamic range of die 
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Figure 13, HERITAGE data lor ilie LMC with the ^/^/--SAGE-LMC MIPS 70 |;m image (Meixner et al, 2006) for comparison. Together these data sets sample ihe 
spectral energy distribution of ISM dust emission. The gray-scale displays are all in aresinh with the following ranges: MIPS 70 pm, 0-400 MJy sr 1 : PACS 100 /*m T 
0—400 MJy s7 PACS 160 //m. 5-500 MJy sr SPIRE 250 /my 1-300 MJy sr J ; SPIRT 350 pin. 0-1 70 MJy sr -1 ; SPIRE 500pm, 0-70 MJy sr L 


Herschel images is larger than that of Ihe II i because the dust 
emission is very sensitive to temperature, ihe bright, hotter 
dust regions in die Herschel images appear morphologically in 
lire same places as the brightest Hty emission in the SHASSA 
images. For example* the brightest Hn regions in the SMC* 
N66* and in the LMC* 30 Doradus, appear as bright regions 
of emission in the PACS 100 pan and SPIRE 250 /an bands. 
Massive 013 stars that ionize gas are also the most effective at 
heating the dust, explaining this morphological correspondence. 
Skibba et al. (2012) presents a detailed comparison of the stellar 
and dust luminosity. 

The CO ./ = 1-0 emission bears little overall resemblance to 
the Herschel bands because most of the ISM gas in die LMC and 
SMC is neutral atomic gas. However, if we look at die structures 
delected by the CO J — 1-0 emission and search for diem in 
the Herschel images, we do find enhanced emission. The ISM 
gas regions detected in CO are associated with giant molecular 
cloud complexes (e g.* Kawamura et al. 2009) and dins with 
larger column densities of ISM gas and emission* explaining 
this correlation (cf,* Bernard el al. 2008), 

Therefore the Herschel emission is affected by both die 
temperature and column mass density of the ISM dust. We have 
presented analyses of the ISM dust mass using the LMC SDP 
strip in Meixner et al. (20 1 0), Gordon et al. (20 1 0), and Galliano 
et al. (201 1 ), Analysis of ISM dust mass and temperature based 
on the Herschel emission from die whole LMC and SMC 
will be covered in a forthcoming paper (K. D. Gordon et al. 
in preparation). Analyses of the dust-to-gas mass ratio in die 


ISM on tlie LMC SDP strip were presented by Meixner et al. 
(2010) and Roman- Duval et al. (20 1 0). Complete analysis of the 
dust-to-gas mass ratio across the whole LMC and SMC by 
comparing the Herschel data widi gas tracers will be presented 
by J. Roman- Duval cl al. (in preparation), 

6. GLOBAL SPECTRAL ENERGY DISTRIBUTIONS 
OF THE LMC AND SMC 

We calculated die total galaxy dux in the LMC and SMC by 
summing over die duxes of die background- subtracted images. 
We then adopted a signal -to-noise (S/N) cut for the pixels 
incorporated into the dux summation* w hich primarily improved 
the PACS 100 /on ilux for the SMC. Table 4 lists the total 
galaxy fluxes for die LMC and SMC for all live bands of die 
HERITAGE survey. We compared these 11 LRU AGE fluxes, 
shown as red triangles, lo prior measurements, w hich are shown 
in blue circles* in the SED in Figure 1 9. The duxes for both 
galaxies measured in prior work are taken from the recent 
compilation of measurements from UV to radio wavelengths 
including measurements from a combination of ground-based 
telescopes (e.g., at optical and liear-IR wavelengths) and space- 
based missions such as IR IS, die Infrared Space Observatory, 
and Spitzer (Israel et al. 2010). f or the SMC, we also include 
the SAGE-SMC MIPS measurements (Gordon et al. 20 1 1 ). 
In order to investigate the submillimeter excess emission in 
the Magellanic Clouds. Israel et al. (2010) measured the total 
fluxes of the LMC and SMC trotn two all sky CMB satellites 
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Right Ascension (J2QG0) 

Figure 14. LMC HERITAGE data. Red corresponds to SPIRE 250 /jriv <>reen to PACS 160 ^ny and blue to PACS 100 jt m. 


Ta b le 4 

Total Galaxy Fluxes in SLD 


Parameter 

PACS X (//m) 


SPIRE X (jum) 



100 

160 

250 

350 

500 

LMC flux (Jy) 
SMC flux (Jy) 

1.6 ± 0.2 x 10* 
1.9 ±0.3 x 10 j| 

2.7 ± 0.4 x 10* 
3.0 ± 0.4 x 10 1 

1.4 ±0.1 x 10* 
1.0 ±0.1 x 10 4 

6.9 ±0.7 x 10 1 

5.9 ±0.6 x 10* 

3.0 ± 0.3 x 10 4 
2.9 ± 0.3 x 10* 


COBH (Roggess et al, 1992) and the Wilkinson Microwave 
A fiksotropy Probe (WMAP\ Bennett et al, 2003a, 2003b), For 
ihe CUBE sa tell lie. they used data Ifom the D1RJ3E instrument 
(Silverberg et al. 1993) and die FIRAS instrument (Fixsen et al. 
1994; Wright et al. 1991). On this topic of the submillimeter 
excess, the Planck Collaboration et al. (201 1 ) reported the flux 
measurements tor the LMC and SMC at Planck wavelengths 
for the total flux, the CMR subtracted flux and the CMR+MW 
foreground subtracted flux. By analyzing these different fluxes. 


they claim that the submi dimeter excess in these galaxies may 
well have a significant contribution from CMR fluctuation 
emission and in the case of flic LMC it may explain all of 
flic excess emission at wavelengths much longer than probed 
by HerscheL In HERITAGE we probe a different type of 
excess emission from the ISM dust (e g., Gordon et al. 2010), 
Our preliminary analysis of the submillimeter dust emission 
measured by HERITAGE suggests the 500 gm excess emission 
regions are correlated with LMC ISM structures so likely 
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Figure 1 A HF RTTAGF d;il.i Ii?r Ihc SMC wilh the SpifzerS AGE -SMC MIPS 70 pm inujiy (Gordon ct a\. 201 1 ) for comparison. Tmpedicr these dal j mjI.s sampl* 
the spectral eneigy distribution of ISM dust emission. The gray-scalo displays are all in arcsinh with the tollraving ranges: MIPS 70 pm, 0 100 MJv sr -1 : PACS 
100 pm 5-250 ALU sr ' ; PACS 160 pm 5-200 MJy sr 5 ; SPIRE 250pm 0^1 00 MJy *r 3 j SPIRE 350 pm 0-SQ MJy sr 1 ; SPIRE >00 pm 0-30 MJy sr j . 



Right Ascension (J2000) 

Tiguri* 16. SMC HERITAGE data. Red corresponds to K1‘LRL' 250 pm, green to PACS 160 pm, ami blue to PAC S 100 pm. 
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Figure 17. LMC survey data. From lop to bottom, left to right: the Spitz? r JR AC" 1 3.6 image from the SAGE-LMC survey (Meixner et al. 2006) reveals the old stellar 
population, which is ejecting dust into the ISM. The Herschel PACS 160 and XP1RF 250 ;im image from the I HERITAGE survey reveal the coolest dust emission, 
which is dominated by the huger dust grains. The [Jo emission from SUAKSA (Gaustad el al. 2001) shows the distribution of the diffuse ionized gas. punctuated by 
11 ii regions. The II i 21 cm emission (Kim et al. 2003) show the location of the atomic neutral hydrogen. The CO./ — 1 -0 emission survey by NANTFN (Fukui el al. 
2008) traces (lie highest column density of molecular gas. The gray scale displays are alt in ancsinh with the following ranges: 1RAC 3.6 p ni 0 5 MJy sr 1 . FACS 
160 p m 5 500 MJy sr” 1 ; SPIRE 250 pm 1 300 MJy sr” 1 ; Ho 0 8000 ; H r 21 cm. 4 70 x 10 20 cm” 2 ; CO 0, 0 1 5 K km s" 1 . 


originate from dust emission in the galaxies (eg., Galliano et al. 
2011; Verdugo et al, 2013). Further analysis of this 500 ptm 
excess emission across the LMC and SMC will be discussed in 
a future paper (K. D. Gordon el al. in preparation). 

In f igure 1 9, we plol Ihe total LMC and SMC lluxes without 
any CMB subtraction measured by the Planck Collaboration 
et al, (20 1 1 ) in order to be consistent with onr HERITAGE 
fluxes and the prior work. We calculate die Planck total fluxes by 
using the reported average surface brightness times the circular 
area defined by a radius listed in the Table 2 caption from 
Planck Collaboration ct al (2011). The Planck data, which 
arc shown in green squares in Figure 1 9, appear in reasonable 
agreement with the prior work that includes the WMAP and 
CY95E/D1RBE/FIRAS results. In comparison to this complete 
SFD of the LMC and SMC, the HERITAGE SPIRE photometry 
data agrees within the respective errors of the data sets The 
HERITAGE PACS 100 /tun data appears to also agree with prior 
work widiin the enorbars. t he HERITAGE PACS 160 fan flux 
point appears consistent with prior measurements for the LMC 
but appears to slightly exceed the known photometry Lor die 
SMC (Gordon et al. 20 1 1 ; Leroy et al. 2007). 

7. HERITAGE POINT SOURCE EXTRACTIONS 

The variable background emission present in the HERITAGE 
mosaics presents a challenging environment in which to attempt 


point-source photometry . After experimenting with several 
source extraction packages, we chose the PSF-fitting software 
Starfmder (Dio laid et al. 2000; see Section 7.3). Aperture 
photometry is not feasible for automated production of catalogs 
because the diffuse emission is too complex for automated 
annulus selections to calculate background levels. Stall inder 
uses a smoothing algorithm to estimate die local background 
and dien iterates on the background as sources are found and 
extracted (Berlin & Arnolds 1 996). 

The HERITAGE catalogs were produced through an itera- 
tive process. T he first step was to find and extract sources with 
Starfmder, using PS I s provided by die Herschel Science Cen- 
ter. 3 * The detection threshold was set at 5a above die noise 
level. The correlation parameter, a measure of goodness of lit 
(where I is perfect), was used to remove sources with values 
below 0.75. Source lists were made from not only the combined- 
epochs image but also each individual epocli image (two epochs 
for SPIRE LMC and SMC data and PACS LMC data, and diree 
epochs for the PACS SMC data). These additional single epoch 
sourcelists were used to improve the reliability of the filial cat- 
alog of sources, which contains only sources detected in the 
combined-epoch images. For each of the live wavelengths of 
data (PACS 100/un and 160/un and SPIRE 250 350 /an, 


b [ Ip: h erse he E .esae .csa . m I twi ki bi n vi e vv Publ i c 1 PacsCa I i bra E i cm Vv 'civ a ml 
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Figure 18. Same type of survey data as in Figure 1 7 but from different sources, for the SMC. The Spifzer IRAC 3.6 ft m image from the SAGE-SMC survey (Gordon 
et al. 201 1 ) reveals the old stellar population. The JJerscheJ PACS 160/nn and SPJRF 250 /on images from ihe HERITAGE survey reveal the coolest dust emission. 
The Ilu emission from SHASSA (Gauslad et al T 2001) shows the distribution of the diffuse ionized gas and I In regions. Hie Hi 21 cm emission (Slanimitovic 
et al. 1999; Muller et af 2003a) show the location of the atomic neutral hydrogen. The CO J — 1-0 emission survey by NANTEN (Mizuno el aL 2001) traces ihe 
molecular gas. T he gray scale displays are all in arcstnh with Ihe following ranges; IRAC 3.6/im. 0 5 MJy sr _l ; FACS 160 /,im. 2 200 MJy sr -1 : SPIRE 250 //m. 
0 100 MJy sr" 1 ; He* 10 5000 : Hi 21 cm, 5 100 x 10~ u em“ : ; CO J ^ 10, 0-2.2 K km s _1 . 


and 500 /4rti) the individual epoch source lists and combined- 
epoch image source list were matched using a simple nearest 
match criteria of 0.7 x FWHM. llie subset of matched sources 
for each band was then identified tor enhanced study lo de- 
termine source reliability and photometric accuracy For these 
subsets we fit a 2D Gaussian to each source. We then applied 
an additional cut to our subset of sources, selecting only the 
sources that produced a measured FWHM in both x and v di- 
rections without extreme values (between l /3 and 3 times the 
Hers die l FWHM value). A histogram of those FWHM values 
for each wavelength band was then used, as well as ihe mam 
and median values lo detenu ine an optimal FWHM for the 
extraction PSF. The Herschel Science Center PSFs were then 
broadened with a Gaussian of appropriate width to produce a 
new PSF with the desired FWHM (see Figure 20). Many factors 
may have contributed to make the HERITAGE PSFs more broad 
than the Herschel PSFs. The HERITAGE mosaics were created 
from lengthy scanning observations, in which pointing uncer- 
tainties. changing position angle, and multiple observations can 
result m a broadened PSF. Unfortunately, the degree of broad- 
ening was variable across images. The PSF from the combined 
image of all the epochs was more broad and circular than those 
from the individual epochs. 

Using these broadened PSFs, Starfinder was re-run on the 
HERITAGE images. An iterative process was employed, llie 
lirst pass of Slariinder extracted llie vast majority of sources. 


Then the residual image (the image after Starfuider removes the 
found sources) was resubmitted lo Starfuider a second time to 
find additional sources that may have been lost in the wings 
of the PSF of bright sources. Additionally, ailer each Slariinder 
iteration, we used a tweakmg algorithm on the residual images to 
correct for over-extracted sources in high /complex background 
areas. This algorithm is described in more detail in die Spine r 
GLIMPSE processing documents. 39 Ibis additional tweaking 
step limits die occurrence of multiple source extraction from an 
overly broad source with a PSF smaller than die FWTIM of the 
broad source. 

7./. Flux Uncertainties 

Starfuider produces an uncertainty associated with each 
extracted source, which is determined by the mode in which 
Starfinder is ran. For the SPIRE data we used the associated 
err or maps produced from die image mosaicking process. For the 
PACS images we allowed Starfinder to detenu ine the uncertainty 
from the nns in the background part of the image, i he PACS 
images contain a degree of low level striping that adds more 
uncertainty to die lower background level sources. I lowever. the 
uncertainty dial Starfinder detenu i nes is only one component of 
several that contribute to the overall uncertainty of the duxes in 
die HERITAGE catalogs These additional components include: 

39 h I Ip: WWW. as I ro.wi sc. cdu gti m psc gl i mpsc photonic! iy v \ .(1. pdf 
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Figure 19. Spectral cncigy distributions of the whole LMC and SMC cover the UV Io radio wavelengths. Prior work summarized by Israel el al. (201 Q) includes 
data from ground and space based missions such as IRAS. DIRBE, TopIIat and WMAP. for the SMC* we also include prior work (blue squares) horn Gordon et al. 
(201 1 ). The HERITAGE photometry measured in the PACS and SPIRE bands are shown in red. The Planck measurements are shown in green (Planck Collaboration 
et al. 201 IX 

(A color version of this figure is available in the online journal.) 


determination of the background* the accuracy of the PSF used 
in the extraction* and the degree to which die extracted source 
is a true point source. 

lliese other components often dominate in die l IERITAGE 
data. We attempt lo quantify these components and give more 
generalized uncertainties to the HERITAGE source lists by ana- 
lyzing die flux measurements of many sources observed in two 
different epochs. Ideally, observing a single source many times 
will provide an accurate estimate of die uncertainties. We ob- 
served each source only two times, typically, but we can compare 
hundreds of different sources of similar brightness and back- 
ground levels. Thus, we matched sources from the combined- 
epoch images to die epoch 1 and 2 images, grouped sources of 
similar brightness and background levels, and compared their 
fluxes using the formalization found in die SAGE-d ,MC Data 
Delivery document (Section 3.5) for dieir parameter f 1u : 


K+*l) 1/2 ’ 


4 ' } See hi i p; sage, slsei .edit S AGE -LMC -del iv2 .pdf. 


( 1 ) 


where F \ is the epoch 1 dux of a source, F? is the epoch 2 flux of 
die same source, and 04 and are the uncertainty estimates of 
die epoch 1 and epoch 2 duxes. The uncertainties are calculated 
as follows, using 04 as die example; 

(T, = (■'f+Sl 2 / + 4) 1 ' 

(2a) 

s\ = unci * uca 

(2b) 

$1 f ~ * /cor 

(2c) 

i 

* 

t 

!! 

(2d) 


where uncj is die output uncertainty from Starfinder. We include 
additional uncertainties s\ f and that depend on die 11 ux 
and sky values (with, coefficients / CjW and x car respectively). 
We also include an additional coefficient n cof on s\ to account 
for additional low 11 ux level uncertainty that the Starfinder 
uncertainty underestimates. Values of H cor , / cor and Scor are 
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PACS100 Combined Epoch image 



PACS160 Combined Epoch image 



FWHM (orcsecs) 


FWHM (orcsecs) 




FWHM (orcsecs) 


Figure 2®. Plotted are histograms of the FWHM (average of the .v and y directions) of HERITAGE catalog sources. The black histograms show all sources that could 
be lit with a Gaussian in both .v and y directions. Sources with the Lowest sky values (40% of loial) are plotted in red. This shows that sources on higher background 
sky areas tend to produce higher FWHM values. This can be attributed to extended emission of the source, or to the more difficult task of separating the point source 
horn the background emission. The black vertical line is the FWHM (2D average) of the initial H ^rsc/ire/suppl ied PSFs, and the red vertical line is the FWHM of the 
broadened PSFs that were used to produce the HERITAGE catalogs. The SPIRE 500 pin PSF values suffer the most from background sky contamination. 

(A color version of I his figure is available in the online journal.) 


empirically determined to optimize the resultant uncertainties 
so dial lire distribution of flux differences V between epoch 1 
and epoch 2* for sources in similar groupings* forms a Gaussian 
distribution with <r — 1 (see Figure 21). Note that the center 
position of 0 in the histograms shows that there are no systematic 


offsets in die duxes between epochs 1 and 2. The values of n coo 
/cor and Scar vary with wavelength and each galaxy (LMC and 
SMC); for example the LMC SPIRR 250 jim single epoch values 
are u ^ — 1.4, = 0.06, ,v COT = 035. Table 5 lists the final 

values for both the LMC and SMC datasets. 
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Figure 21. Each panel shows a subset of sources from the SPIRE 250 /rm catalog plotted as histograms of the variability parameter F (Equation (1)). Eaeli pane! 
has an equal number of sources but has distinct ranges of sky and source Duxes, which are labeled on the left and right panels amt are intermediate in values for (he 
middle panels, lire Dux ranges increase from left to right, and the sky values inerease from lop So bottom. The red. thin eurve shows the expected histogram width if 
(he distribution of i s has a — 1. A histogram broader than (he red curve implies unceilainties are loo small compared to the difference between Dux values (assuming 
minimal number of variable sources). The green, thickest line histograms show the values based on Slarlinder unceilainties only. The width of the green histogram 
increases for both increased sky and Dux values indicating larger differences, and therefore larger uncertainties, in the epoch 1 and 2 (lux measurements. The blue, 
thinner line histograms show the values based on the adjusted source uncertainties as described in the text. These blue histograms are much closer in value to the 
distribution of V% such that tj — 1. In our analysis, we maximized the number of grids for better differentiation when jrossible. but here only presenl the minimum 
3x3 grid, where each grid contains 1620 sources. 

(A color version of this figure is available in the online journal.) 


Additional uncertainty exists in the FACS source extractions 
that is not accounted lor in die variability parameter analysis 
above which results from the l/f noise in masked regions, 
lire 1 / f noise in single epoch images causes striping in one 
direction, while combined-epoch images have striping in two 
perpendicular directions, causing a cross-hatching pattern. Tire 
noise is higher in the cross-hatched regions of the combined 
image with stripe intersections coadding to increase the noise. 
Wo isolated die hatching in order lo measure the magnitude of 
diis added uncertainty by subtracting die Starimder-produced 
sky image from die residual image, thus removing die ISM 
structure from die residual image and leaving only the cross- 


hatching pattern and other instrumental noise. Point sources 
and ISM emission have been removed, so the flux within a point 
source-sized aperture should be zero. We calculated the flux 
within 10.000 apertures (aperture radius = i?so% ; background 
annulus inner radius = 1.2 x Rm%- annulus w idth = 2 pixels, 
where is the radius which contains 80% of the psfs 
power) in both masked and unmasked regions of both galaxies' 
FACS images and measured the scatter about a flux of zero. 
The masked regions, which suffer strongly from striping, have 
an enhanced scatter (Figures 22 and 23), indicating a larger 
II ux uncertainly for sources located in masked regions relative 
to those in unmasked regions. To account for this additional 
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PACS 100/im PACS 160^im 



F ioo/im [mjy] [mjy] 

Figure 22. Histograms showing the distribution of the aperture photometry of 10,000 apertures on the LMC PACS residual minus sky image. T!ie histograms are 
normalized to peak at /V = 1. PACS 100 /jm is shown to the left and PACS 160 jum to the tight. The sealier of the 13 ux (quantified with the standard deviation, cr) is 
indicated in each panel for apertures in unmasked {Sop row) and masked (bottom row) regions. 


Ta blr 5 

Uncertainty Correction Factors for the Combined-epoch linages 


A fern) 


feat 

s cat 

LMC 

100 

3.00 

0.065 

0,25 

160 

400 

0.060 

025 

250 

1.40 

0.060 

0.35 

350 

1.60 

0.050 

120 

500 

T50 

0.070 

2,50 

SMC 

100 

3.00 

0.050 

0.40 

160 

3.00 

0.065 

0.55 

250 

1.25 

0.030 

0.50 

350 

1.20 

0.050 

1.20 

500 

L45 

0,070 

2.40 


uncertainty from the 1 // noise, we have linearly added an 
uncertainty of (rr masked - rr imm ^\:ed) to the uncertainty determined 
from tlie variability parameter, where amukad is the standard 
deviation of the scatter about a Jlux oi zero in masked regions, 
and ^unmasked is the same in unmasked regions. In the LMC, 
we have increased die PACS 100 /an and 160 /an uncertainties 
by 6 mJy and 17 mJy, respectively. For the SMC. die same 
uncertainties are increased by 6 mJy and ! 1 m Jy. Recall that die 


difference in sensitivity between masked and unmasked regions 
is greater for PACS 160 /an than for PACS 100 /an (Table 3), 
consistent with the larger value of (Gm&ked ~ ^unmasked) for PACS 
1 60 /an than PACS 1 00 /an. 

The scatter about zero flux from this same "source- less 7 " 
aperture analysis in SPIRE and unstriped PACS regions is 
a measure of the minimum uncertainty that results from the 
instrumental noise of the image. This is essentially the minimum 
value of s (and therefore o ) from the variability parameter, as 
a = ,s when A and sky are both 0. The minimum uncertainty 
derived from the aperture analysis and that which results 
from the multiplication of unc and show good agreement 
(within 1-3 mJy). Perfect agreement is not expected, as one 
is determined from an aperture photometry analysis while the 
odier from Starfinder, but their close correspondence confirms 
reasonable values of u txjr . 

in Figure 24 we plot the uncertainty of the lluxes as a function 
of flux for so mces extracted from die PACS 100 // m and SPIRT 
250 /an combined-epoch images. Lines of constant S/N are 
drawn to guide die eye. Note that an upper limit of the S/N is 
set by /cor at higher fluxes and u cor at low fluxes. 

These uncertainty values provide an accurate assessment of 
the uncertainties of the lluxes in die HERITAGE catalogs. For 
detailed examination of particular regions, manual photometry 
on individual sources may improve the flux estimates, by 
allowing the user to vary the PSF, apply non-circular aperture 
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PACS i00/im PACS 160/mi 



Fioo„m f m JyJ F ieo^m [mjy] 

Figure 1 23. Same as Figure 22 . but for the SMC. 


photometry, or heller determine the local background for each 
source. In the case of PACS, when possible, locally de striping 
the image before attempting photometry may improve both die 
flux and uncertainty . 

7.2. Simulated Source Tests 

A battery of simulated data tests was run on both die PACS and 
die SPIRE LMC data sets to assess completeness and photomet- 
ric fidelity. A 5 deg 2 area centered at 85:7 6 9 8 (R.A., dec!) 

was selected based on its wide range of background variability, 
lire region was populated with a selection of 800 simulated 
sources. Both die combined-epochs image as well as the single 
epoch images were tested. The 800 sources consisted of eight 
flux bins each with 1 00 sources. The flux bins were 20,50, 100* 
200, 500, 1 000, 2000 and 5000 mJy. Source locations were cho- 
sen to avoid locations within 2 xFWHM where known sources 
were extracted. We instituted this source location restriction 
so that we could assess the completeness as a function of back- 
ground and not confusion from crowding. In diese datasets, con- 
fusion widi other sources is low compared to confusion with die 
high, complex background. For SPIRE, sources were inserted 
into the timeline datasets and subsequently those timeline data 
were re- mosaicked into map images from which the photometry 
was done. This was done so dial we could simultaneously deter- 
mine the value of the pixel Nation correction (SPIRE Photome- 
try Cookbook 51 ). Pixelization correction values of 0.95 1 , 0.93 1 
and 0.950 were determined for SPIRE 250, 350, and 500 /an 


11 litlp: hcrschd.csac.csa.int (wiki pub/Public^SpireCali bra lion Web 
SPUUiPliolomcIrv 'Cookbook jiprZOl 1 _3.pdf 


data respectively. Willi die pixelization correction applied to 
die SPIRE results the photometry of simulated sources rqiro- 
duced true fluxes well No pixelization correction is required 
for PACS data, so sources were inserted directly onto die map 
images. 

fhe simulated source tests found a tight agreement between 
true and extracted source fluxes. The ability of the photomet- 
ric process to find and extract all simulated sources (the mea- 
surement of completeness) is a strong function of the back- 
ground flux and complexity. For example, die left-hand panel of 
Figure 25 shows die completeness curve for sources on a low 
background, while die right-hand panel is die completeness 
curve for regions we define as high background (see Table 6 
for the del hi i bon of high- and low-backgrounds), fhe overall 
completeness at a 90% level, as determined by these simu- 
lated data tests, is tabulated in fable 6. lire completeness is 
clearly a function of background, and given the highly variable 
nature of die background, die completeness has substantial vari- 
ability even within the high-background regions. These curves 
therefore represent the average completeness within die region, 
widi die lute completeness of a particular region depending on 
the brightness and complexity of die background immediately 
surrounding the source. Global LMC and SMC average 90% 
completeness limits are calculated for each wavelength taking 
into account the 111 ling- fraction of the galaxy with high and low 
surface brightness fluxes (Table 6) These average values under- 
estimate the catalogs' completeness levels because while most 
of the image is low background, most sources are located in 
high background. Nevertheless, they provide rough guidelines 
for die limitations of the catalogs. 
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Figure 24. PACS 100 /on (lop row) and SPIRE 250 jum (bottom row) catalog source llux uncertainly ;is a function of H ux for ihc LMC (left). and SMC (right). Lines 
of cons la nl S/N are drawn to guide Ihe eye. 


Low Background High Background 



Figure 25, Completeness curves for SPIRE 250 pm created from the false source injection test, which show the percentage of false sources recovered as a function 
of flux, ihe left and right panels are for sources within a low (^10 MJy sr _1 ) and high (> 10 MJy sr” 1 ) background, respectively. The black curve is for all sources 
recovered from Ihe combined image (i.e., Ihe Pull List), while the gray curve is Full List sources delected in both individual epoch images. The vertical lines mark ihe 
Full List J?0% completeness limit at each background level. 


7. 3. Soane Extraction Assessment and Applicability (Mo I man et al 20 1 ] ), and IRAFs Apphot (Tody ! 993 ). Each of 

these packages extracts tire flux and detenu ines die background 
in a slightly different manner; Starfitider uses PSF fitting* 
CuTEx uses a 219 Gaussian fitting routine, and 3 Extractor 
and Apphot arc aperture photometry routines. Wc ran each 


To verify that Star! inder w as producing reasonable results, wc 
performed rough extractions using alternate source extraction 
packages including SExlraclor ( Berlin & Amouts 1 996), C uTEx 
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Table 6 

Parameters of Extracted Point-source Catalogs 


Parameter PACS A. (/im) SPIRE a (jum) 

100 160 250 350 500 


LMC 


Total number in extracted Fist 
Number in catalog 
FWHM of extracted source (") 

86 86 
4164 
8.6 

27136 

9324 

12.6 

48082 

25445 

18.3 

38482 

22082 

26.7 

12224 

7355 

40.5 

SMC 

Total number in extracted list 

1353 

5810 

11923 

11669 

2335 

Number in catalog 

898 

1590 

5465 

5313 

1069 

FWHM of extracted source (") 

8.8 

12.6 

18.3 

26.7 

39.7 

90% completeness limits from simulated source tests 

Low background (MJy sr l ) 

^15 

<25 

<10 

^5 

<2.5 

\ .ow background 90% completeness limit (mJy) 

450 

150 

50 

50 

90 

l ligh background (MJy sr 1 ) 

>15 

>25 

>10 

>5 

>2J 

High background 90% completeness limit (mJy) 

450 

400 

300 

400 

400 

LMC average 3 90% completeness limit (mJy) 

450 

190 

90 

no 

140 

SMC average 3 90% completeness limit (mJy) 

450 

160 

60 

60 

100 


Note. a Hie average is computed as Ctotal ” Clow * /low + <^high * /hi£h? where Cis Ihe 90% completeness limit, /is the filling file lor 
across the galaxy, and the subscripts tow and high reference the corresponding background level. 


package on a subregion of the LMC that contained both high- 
and low-background regions to gauge the consistency between 
the methods, and found each package to return source fluxes 
that generally agreed to within a lac tor of -'1.5, with the 
largest discrepancies for sources in lire brightest, most complex 
backgrounds. Sources in more isolated regions of the SPIRE 
250 /t in map generally agree to within 20%. The tests confirm 
that Starflnder returns reasonable results, and affirms that one 
of the primary sources of flux uncertainty is the variable 
background. 

As described in Section 7.2, simulated source tests revealed 
a tight agreement between the Starfl nder-extrac ted and true 
fluxes, t his is perhaps unsurprising as the inserted sources 7 
PSF matched that used during extraction by Starflnder. For 
real sources, this is not necessarily the case; indeed, we find 
some variation in PSF width for the HERITAGE sources (see 
above). The Starflnder post-extraction residual images (i.e. : the 
original image minus the PSF fits) show some sources with 
unextracted, residual flux at the position of a source. Hiesc 
sources have extended emission beyond the PSF such that 
while the bright point-like center of the source is removed, 
the lower level, more extended flux is not. The inclusion of this 
residual flux in the source's flux depends on die interpretation 
and science goals of the user. Many of the sources identified 
in the HERITAGE images are embedded in large-scale dusty 
structures, and the sources themselves may not be discrete 
compact structures, but rather localized bright, dusty density 
enhancements. Such a source may not be strictly point-like, and 
aperture photometry may be better suited for such sources, but 
decisions about using PSF-fitled or aperture photometry must 
be made on a source-by- source basis. Moreover (lie varying 
PSF width may require a variable aperture correction factor. No 
single photometric method will fill all desired photometric roles, 
however we found the Starlinder results produce good results 
for true point sources and lower bounds for sources within areas 
of high complexities. To designate good point sources we fit 
a 2D Gaussian at each source location. To isolate the source 
and minimize any background effects on the 2D Gaussian fit 
we used the residual image (sources removed except the desired 
source) with the background image removed. We tabulate the 


Ta bic 7 

Flux Hag Statistics: Units arc Percentage of Sources 
in ihe Catalogs with Each Flux Hag Value 


U/im) 

Hag = 1 

I-lag = 2 

Cl 

III 

Hag = 4 

LMC 

100 

44 

9 

7 

39 

160 

31 

4 

12 

54 

250 

35 

16 

34 

13 

350 

29 

20 

26 

24 

500 

21 

23 

39 

16 

SMC 

100 

73 

5 

7 

15 

160 

66 

2 

7 

24 

250 

74 

4 

17 

4 

350 

65 

10 

20 

5 

500 

54 

12 

26 

8 


results as a FWHM value in our catalog (averaging the _r and v 
FW1 1M values), sources where the Gaussian fit haled (in either 
direction, or both) are designated with the value of 0. 

These FWI1M values determine a photometric quality flag, 
flux flag, 7 in our catalog indicating the degree to which die 
source is point-like and thus the reliability and applicability of 
die flux measurement returned by' Starflnder. The flag has a value 
ofl , 2, 3, or 4, and Table 7 documents how many sources in each 
galaxy and waveband are assigned each flag. A flux flag of 1 is 
given to sources that have a measured FWHM consistent with 
the FWHM of Are PSF (FWllM B(HirCL . < 1 . 1 x FWUMpsp). T he 
flux for flux flag = 1 sources is therefore a robust accounting of 
the sources’ brightnesses. 

A flux flag of 2 is assigned to sources for which the FWHM 
detenu i nation failed, but die source has a S/N of greater than 5. 
Tire great majority of flux flag = 2 sources are faint, causing the 
FWHM measurement to fail from a poor Gaussian fit. Bright 
flux flag = 2 sources tend to tail FWHM measurement due to 
their location within areas of complex background. For SPIRE, 
this constitutes a small fraction of sources: 0.6% of 250, 350, 
and 500 /nn sources are brighter than 0.5 Jy and have a flux 
flag of 2. Because the PACS sensitivity is poorer fewer faint 
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sources have a S/N > 5 (one of the llux llag=2 requirements), 
and thus there is a higher fraction of bright llux Hag = 2 sources; 
^4% of PACS sources are brighter Ilian 0.5 Jy and have a flux 
flag of 2. 

In Figure 26* we present small image cutouts of the area 
around typical sources identified in the SPIRE 250 /um image. 
The sources are organized by flux flag and flux. Note that 
faint flux Hag = 2 sources are located in relative isolation 
and that the residual images show little to no signs of under 
or over extraction. Conversely, bright sources are located in 
regions of more complex background and often show signs of 
residual flux surrounding lire source. In these cases, Starfinder 
lias extracted the central point-like core of the source, leaving the 
more extended surrounding envelope* and whether the extended 
emission should be included in the source flux depends on tire 
science case. We leave it to the user to determine if the extended 
emission should be included in the source's photometry, in 
which case the user may wish to increase the uncertainty of 
the Starfinder extracted llux or treat it as a lower limit. Again, 
this constitutes a very small fraction of the Catalog, so most flux 
flag = 2 source fluxes can be treated as good measurements. 

Flux flag of 3 denotes sources that have a measured FWHM 
greater than that of the PSF (FWHMgoqree > LI x FWHMpsf) 
and a S/N > 5. These sources are extended beyond the PSF, 
and thus Starfinder has extracted the point- like central core of 
flic source (see Figure 26). In some eases tire user may wish to 
include this extended flux, in which case aperture photometry 
may be more useful and tire Starfinder flux can be considered 
a lower limit Sources with a S/N < 5 are marked with a 
flux flag of 4, and have uncertainties that are dominated by 
the brightness and complexity of the background and for some 
PACS sources l// noise; thus the extracted flux is highly 
uncertain. Note in F igure 26 the highly complex nature of these 
sources' environments. Because the sources in the catalogs we 
provide the community have passed fairly rigorous tests to verify 
lire sources are real (see below), flux flag — 4 sources should 
be considered tme sources with highly uncertain or mi known 
fluxes The fluxes provided in die catalogs for flux flag = 4 
should be treated with caution, and a user may wish to do a 
more detailed source-by-source extraction to obtain a higher 
fidelity measurement. 

Note from fable 7, that there is a larger percentage of flux 
Jlag — 4 PACS sources than SPIRE sources, lliis is a result of 
the poorer PACS sensitivity, particularly in the masked, striped 
regions where a large fraction of the sources reside. There is 
also a higher percentage of flux flag — 4 sources in the LMC 
than die SMC due to the LMC’s higher sky levels. Finally, note 
that the number of flux flag = 1 sources generally decreases 
with increasing wavelength while the number of flux flag = 2 
or 3 sources increases. As the observed wavelength increases, 
the emission surround big point sources originates from cooler, 
more extended dust, thus broadening the FWItM and relegating 
more sources to flux flag = 2 or 3 status. 

7.4. Assessing she Reliability* of Sources 

Table 6 documents the number of sources identified in each 
combined-epoch image by Starfinder. These sources represent 
the “Full List” of point sources, every source identified in die 
combined-epoch images. We provide tor use by the community 
a subsample of die l ull List that has been vetted for reliabil- 
ity, which w r e call the “Catalog.” The Catalogs contain only 
combined-epoch sources that pass die reliability tests detailed 
below The Full List may contain sources — particularly at low 


flux levels — dial are false detections. Moreover, because of the 
bright, highly variable, and complex background present in our 
images, it is necessary to assess the point- like nature of sources 
to separate tme point sources from marginally point-like fluc- 
tuations in the background, the most reliable point sources 
are those whose point-like emission is quite distinct from the 
background. Table 6 documents the number of vetted Catalog 
sources identified in each combined-epoch image. The HER- 
ITAGE Catalogs are available electronically from the Herschel 
Science Center, and in the Appendix we document the formal 
and column description of the catalogs. 

Our assessment of the reliability of the point sources is a 
two-step process. First we position matched the Full List for 
a particular waveband lo the sources identified in the single 
epoch images of the same waveband. We call tins procedure 
“cross-epoch matching,” A matching radius of 0.7 x FWHM 
is adopted. A source is more reliable if it is identified in both 
the combined image and one or more of the individual epochs. 
False or low S/N sources are less likely to be identified in the 
individual epochs, which are noisier than llie combined image, 
and marginally point- like sources are less likely lo be identified 
as point sources in both the combined epoch and the individual 
epochs. Thus, matching to the individual epoch lists tends to 
identify the brightest point- like sources. Figure 25 displays 
results of die aforementioned simulated source extraction on 
the SPIRL 250 ^m image and compares the percentage of all 
recovered sources detected in the combined image (black line) 
with those detected in the combined image and both individual 
epochs (gray line). The black line shows the percentage of real 
sources recovered during source extraction, and thus represents 
the completeness of the Full List as a function of source flux. 
Note that bright sources are identified in die combined image 
and both individual epochs, while die dimmest sources are only 
detected in die combined image and one or no individual epochs. 

The second step in die source reliability procedure, a source 
position match across ditferent in id- infrared (tnid-IR) and FIR 
wavelengths, is an ellort to identity real sources not detected in 
multiple epochs. We call this procedure “cross-band matching; 
We have position matched the combined-image l ull List from 
each waveband to sources in the oilier Herschel combined-image 
Full Lists and the inid-IR Spitzer SAGE catalogs at 24, 70, and 
1 60 /an. Because of the large variation in angular resolution 
between the wavebands, we only match images with angular 
resolutions that agree within a factor of M.5. For each matching 
pair, a match mg radius of 0.7 x FWHM max is adopted, where 
FWllMjnax is the larger FW r HM of (lie two images. Table 8 
documents the source lists used in the cross-band matching for 
each Uersche! Full List. 

After completion of both of these matching procedures 
(cross-epoch and cross-band matching), we identify the most 
reliable sources Lo be included in the “Catalogs” as follows. 
The reliability is determined by adding n epoch* the number of 
individual epochs a Full Fist source is identified in (n 9p0 oh = 0, 
1, 2, or 3), and /*umd the number of images tor which there is a 
cross-band match (« band = 0, 1,2, or 3). If the sum is 2 or greater, 
the source is reliable and categorized as a “Catalog” source. All 
other sources are less reliable and are included in the Full List, 
which includes every source identified in the combined- epoch 
images. The Full List is not delivered to the Herschel Science 
Center because it is not as reliable as the Catalog; however, it can 
be accessed by contacting the HERITAGE team. The Catalog 
of sources at a particular wavelength therefore contains only 
combined-epoch sources that ( 1 ) are detected in at least two 
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Flux Flag = 4 



Figure 26. Image cutouts showing the combined-epoch SPIRE 250 fim image around a source before (original image; left) and after (residual image; right) source 
extraction. The name of the source is shown above each image pair and the position of the source is marked with a circle. Sources are organized by flux flag and flux 
as shown lo the left. 


29 









The Astronomical Journal, 146:62 (35pp), 2013 September 


Meixner et al. 



Figure 27. Plotted are the distribution of extracted sources as a fund ion of II ux for the LMC (left) and SMC (right). The red cuives are only sources classified 
as Catalog sources, and the black curves are for all extracted sources (Full List sources). Data were binned logarithmically with bin si?.e of 0.1. The average 90% 
completeness level of the catalog (Table 6) is marked by the vertical dashed line. This average 90% completeness level is computed as a weighted average of the high- 
and low -background values and underestimates the real 90% completeness level (see text in Sections 7.2 and 7.4). 

(A color version of this figure is available in the online journal.) 


Table 8 

Cross-band Matching 


Full List Waveband 


fists Matched To 


PACS 100 
PACS 160 
SPIRE 250 
SPIRE 350 
SPIRE 500 


PACS 160 Full List, MIPS 24 Catalog 
PACS 1 00 Full List. SPIRE 250 Full List, MIPS 70 Catalog 
PACS 160 Full List, SPIRE 350 Full Lisl, MIPS 70 Catalog 
SPIRE 250 Full Lisl. SPIRE 500 Full Lisl, MIPS 70 Catalog 
SPIRE 350 Full List, MIPS 160 Catalog 


single epoch images of the same wavelength, (2) are identified 
in at least two oilier wavelengths in combined-epoch Hers cite! 
or Spitzer images, or (3) are delected in one single epoch image 
of the same wavelength and at least one other combined-epoch 
image at a different wavelength. For example* a PACS 1 00 /mi 
source that is detected in two PACS 1 00 /mi single epoch images 


(w epoch = 2) and has a coincident source in the PACS 160 /mi 
combined-epoch image (ft band = l) lias n tpocb + /iband = 3, 
and is a Catalog source. Conversely, a PACS 100 /an source 
not detected in any individual epochs and without a coincident 
PACS 160 /mi or MIPS 24 /mi source is unreliable, a possible 
false detection, and is included only on die Full Fist. Table 6 
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LMC SPIRE 250 Catalog Source Density LMC SPIRE 250 Image 



Figure 2S. Left: ibe l.MC's Kpaiial distribution of the 250 gm catalog sources shown as a source density plot where the grayscale is linear with a range from 0 to 40 
sources per pixel, where ihe pixel is 6 '.6 or 1 00 pc on a side. Right: comparison to the SPIRF 250 gm image shown in gray scale with a square root stretch and a dux 
range of 0 100 MJy sr 1 . 


SMC SPIRE 250 Catalog Source Density SMC SPIRE 250 Image 



Figure 29, LcR: I he SMC's spatial distribution of I lie 250 gm catalog sources shown as a source density plot where I he grayscale is linear wish a range from 0 to 23 
sources per pixel is 5 '8 or 100 pc on a side. Right; comparison to the SPIRF 250 gm image shown in gray scale with a dux range of 0-36 MJy sr 1 . 


documents the number of Catalog sources, that passed all this 
screening, for each combined-epoch image. The Catalog sources 
represent between 27% and 66% of the Full Lisl sources across 
ihe five bands of HERITAGE. 

In Figure 27 we present Hie histogram of the number of LMC 
and SMC sources as a function of flux (in mJy) as well as 
reliability classification (Catalog or Full List). As expected, 
the Catalog In slog ram (in red) euls oil at higher fluxes than 
the f ull Lisl histogram (in black). Using the results of the Lalse 
source extraction procedure described above, we can estimate 
the completeness of the Catalog at low background. Since the 
black line in F igure 25 is the completeness curve for the Full List, 
and the gray curves show the completeness for sources detected 
in both epochs* the Catalog must lie somewhere between die two, 
The Catalog's creation relies on matching across wavebands, 
so the precise shape of die Catalog completeness curve for a 
given population depends on die SLD of dial population. For 
the SPIRE 250 /mi sources, die lesl suggests that the Catalog is 
complete above 100 mJy al low background, with completeness 


quickly dropping below dial level. Recall from above that die 
completeness is highly dependent on the background flux and 
complexity. T lie vertical dashed lines in Figure 27 represent the 
average 90% completeness limits from Table 6 and provide a 
rough guide to the limitation of the catalogs across each galaxy 

8. DUSTY POPULATIONS IN THE LMC AND SMC 

The extracted sources in die PACS and SPIRE bands represent 
some of die dusty inhabitants of die L MC and SMC as well as 
background galaxies. The dux histograms of all die PACS and 
SPIRE bands show that the SPIRE 250 /an band is the most 
sensitive and the PACS 100 /yin is die least sensitive. Thus using 
the SPIRE 250 /an source counts as a guide, we find ^25,000 
sources in the LMC and ^5500 in the SMC. T hese numbers 
are large in comparison to current massive YSO candidate 
lists (M 800 in the LMC and M200 in the SMC) based on 
L Vp?ta?/-only data (see Whitney et al. 2008; Sewilo et al. 2013), 
Figures 28 and 29 display die spatial distribution of the SPIRE 
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250 fim catalog with respect to tlio ISM dust emission at SPIRE 
250 /mi. Hie majority of sources show a strong concentration in 
the bright filaments of dust-emitting ISM; however, there is also 
a background of some uniformly distributed sources. Analysis 
and characterization of these populations will be covered by 
J Seale et al. (in preparation). However, our initial results from 
HERITAGE indicate that those sources associated with the EMC 
are most likely dominated by YSOs (e g., Sewilo et al. 2010) 
and some may be associated with dust clumps (c.g., Kim et al. 
2010). We expect that there are only a few massive and very 
dusty evolved stars (c.g., Boyer et al. 20 10; Clayton et al. 20 1 1 ), 
and one of these is the youngest supernova remnant. SN 1 987A 
(Matsuura et al. 201 1 ). 

9. HIGHER LEVEL HERITAGE DATA PRODUCTS 
FOR THE COMMUNITY 

We have delivered several higher level data products to the 
Herschel Science Center that are available to the astronomical 
community. Combined-epoch images of the EMC and SMC, 
which arc shown in F igures 1 3 and 1 5, arc available in all the 
FACS and SPIRE bands of HERITAGE: 100, 160, 250, 350 
and 500 tim. Hie source extractions are delivered in the form of 
catalogs. Although cross-band and cross-epoch matching was 
used to assess the reliability of the source, the source extractions 
are made available as independent lists for each of the live 
HERITAGE bauds. In the J. Seale el al. (in preparation) paper, 
a cross-band matched catalog will be discussed and published 
for use by the astronomical community. We anticipate that these 
publicly available data products will further the science done 
with llERJ IAGE data, beyond what tire team had envisioned. 

10. SUMMARY 

In this paper, we presented an overview of die Herschel 
HERITAGE project, which surveyed die EMC and SMC in live 
bands including PACS 1 00 and 1 60 fim and SPIRE 250, 350, and 
500 fim wavelengths. We described die observing strategy and 
data processing particularly for PACS, which resulted in maps of 
these galaxies in all live bands. The emission in these Herschel 
bands arises from die ISM dust emission from these galaxies. 
Its spatial distribution is more similar to the ionized atomic and 
molecular tracers than the stellar population. Tile total fluxes 


of the galaxies in these five Herschel bands are reported and 
agree very well with prior total dux measurements, excqit the 
SMC PACS 1 60 fim band slightly exceeds die Spi tzerj MI PS fim 
measurement. We have created source extraction catalogs for the 
LMC and SMC in all five bands. The SPIRE 250 fim catalogs 
are the most sensitive and have —25,000 sources in the LMC and 
—5500 in die SMC. Our Jirst papers, based largely on the EMC 
SDP strip (Meixner et at. 2010), demonstrate the potential ill is 
survey will have for studies in the EMC and SMC. These first 
papers included studies of supernova remnants (Matsuura el al. 
2011; Otsuka el al. 2010), evolved stars (Boyer et al. 2010), 
YSOs (Sewilo et al. 2010; Clayton et al. 2010), II u regions 
(Hony et al. 2010), dust clumps (Kim el al. 2010), dust mass 
and submi dimeter excess (Gordon et al. 2030; Galliano et al. 
20 1 1 ), and dusl-lo-gas mass ratios and related issues (Roman- 
Duval et al. 20 1 0). We anticipate future papers on these same 
topics across the whole of die LMC and SMC. 
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APPENDIX 

MAKEUP OF HERITAGE POINT-SOURCE CATALOGS 

Tables 9-13 document die columns provided in the 
HERITAGE point- source Catalogs and describe the columns' 
contents. The catalog for each waveband and galaxy is pro- 
vided separately; and are available from the Herschel Sci- 
ence Center. Source name convention is an 1 1 character string, 
OOOllPPPPEL, followed by source-spec ill c position designa- 
tion JRR.RRRRSDD.DDDD. OOO denotes die observatory. 


Table 9 

PACS 100 j/m Catalog Contents 


Column 

Name 

Description 

Null 

1 

Source name 

The identifying name of the source. A description of source name construction is found above 


2 

R.A. (J2000) 

Right ascension. J2000 (deg) 


3 

Deck (J2000) 

Declination. J2000 (deg) 


4 

OR A 

Righl ascension uncertainty (aresee) 


5 

dDecl. 

Declination uncertainly (aresee) 


6 

flax 

Star lindcr-exlrae led flux (mJy) 


7 

unc 

Uncertainly of flux (mJy) 


8 

flagFlux 

Hag describing the quality of (lie flax measurement (sec Section 7.3) 


9 

flagS tripe 

Hag notes if source is in masked region and may suffer from striping. ! for masked: 0 for unmasked 


10 

uncSF 

Star lindcr-es Unrated uncertainty of flax (mJy) 


11 

sky 

Starlin del-determined sky brightness (MJy sr 1 ) 


12 

fwhm 

Full-vvidlh-at-half-maximum (aresee) (see Section 7.3) 

0 

13 

flagConf 

Confusion flag (see Section 7), 1 for confused; 0 for not confused 


14 

coir 

Slarfinder correlation parameter (see Section 7) 


15 

epochs 

Number of single-epoch images the source was identified in 


16 

pi 60 

Number of PACS 160/>m Full List sources within 878 


17 

in 24 

Number of MIPS 24 //m Catalog sources within 6:0/6" 2 in the LMOSMC 
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Table 10 

PACS 160 /im Catalog Contents 


Column 

Name 

Description 

Null 

1 

Source name 

The identifying name of the source. A description of source name construction is found above 


2 

R. A. (.12000) 

Right ascension, .T2000 (deg) 


3 

Deck (J2000) 

Declination, J2000 (deg) 


4 

dRA 

Right ascension uncertainly (arc sec) 


5 

dDccl. 

Declination uncertainty (arcsec) 


6 

flux 

Starflnder-extracted flux (mJy) 


7 

unc 

Uncertainly of flux (mJy) 


8 

flaeFlux 

Flag describing the quality of the flux measurement (see Section 7.3) 


9 

flagStripe 

Flag notes if source is in masked region and may suffer from striping. 1 for masked; 0 for unmasked 


10 

irncSF 

SUirfindeT-es lima led uncertainly of flux (mJy) 


11 

sky 

Starfindcr-dctcrniincd sky brightness (MJv sr -1 ) 


12 

fwhm 

FuU-width-at-hall-maximum (arcsec) (see Section 7.3) 

0 

13 

flagConf 

Confusion flag (see Section 7). 1 for confused; 0 for nol confused 


14 

COtT 

Starfinder correlation parameter (see Section 7) 


15 

epochs 

Number of single-epoch images the source was identified in 


16 

pi 00 

Number of PACS 100 /xtn Full T/ist sources within 8'.' 8 


17 

s250 

Number of SPIRE 250 put i Full List sources within 1278 


18 

m70 

Number of MIPS 70 /im Catalog sources within 1276 




Tahle 1 1 




SPIRE 250 /im Catalog Contents 


Column 

Name 

Description 

Null 

1 

Source name 

The identifying name of the source. A description of source name construction is found above 


2 

R.A. (J2000) 

Right ascension, J2000 (deg) 


3 

Deck (J2000) 

Declination, .T2000 (deg) 


4 

dRA 

Right ascension uncertainty (arcsec) 


5 

dDecl. 

Declination uncertainly (arcsec) 


6 

flux 

Statfinder-exti acted flux (tn.Ty) 


7 

unc 

Uncertainty- of flux (mJy) 


8 

flagFlux 

f lag describing the quality of the flux measurement (see Section 7.3) 


9 

flagStripe 

Flag notes if souree is in masked region and may- suffei from striping. 1 for masked; 0 for unmasked 


10 

xmcSF 

Starfinder-estimated uncertainty of Flux (mJv) 


11 

sky 

Slarfinder-delermined sky brightness (MJy sr -1 ) 


12 

fwhm 

Full-width-at-half-niaxitnutn (arcsec) (sec Section 7.3) 

0 

13 

flagConf 

Confusion flag (see Section 7). 1 for confused; 0 for not confused 


14 

corr 

Starfinder correlation parameter (see Section 7) 


15 

epochs 

Number of single-epoch images the source was identified in 


16 

pi 60 

Number of PACS 160 /mi Full List sources within 1278 


17 

s350 

Number of SPIRE 350 /mi Full List sources within 1877 


18 

in 70 

Number of MIPS 70 pat l Catalog sources within 1 278 




Table 12 




SPIRE 350 /mi Catalog Contents 


Column 

Name 

Description 

Null 

1 

Source name 

Ilie identifying name of the source. A description of source name cons Iruc (ion is found above 


2 

R. A. (.T2000) 

Right ascension, J2000 (deg) 


3 

Deck (J2000) 

Declination, J2000 (deg) 


4 

dRA 

Right ascension uncertainly (arcsec) 


5 

dDeck 

Declination uncertainty (arcsee) 


6 

flux 

Starflnder-extracted flux (mJy) 


7 

unc 

Uncertainty of flux (mJy) 


8 

flaeFlux 

Flag describing the quality' of the flux measurement (see Section 7.3) 


9 

flagStripe 

Flag notes if source is in masked region and may suffer from sfriping. 1 for masked; 0 for unmasked 


10 

imcSF 

Slarfinder-es lima led uncertainly of Flux (mJy) 


11 

sky 

Starfinder-detennined sky brightness (MJy sr -1 ) 


12 

fvvhm 

Full-widlh-at-half-maximiun (arcsee) (see Section 7.3) 

0 

13 

flagConf 

Confusion flag (see Section 7). 1 for confused; 0 for not confused 


14 

con- 

Starfinder correlation parameter (see Section 7) 


15 

epochs 

Number of single-epoch images the source was idenlilied in 


16 

s250 

Number of SPIRE 250 pan Full List sources within 1 877 


17 

s500 

Number of SPIRE 500 pan. Full List sources within 2874/2778 in the LMC/SMC 


18 

m70 

Number of MIPS 70 /im Catalog sources within 1877 
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Table 13 

SPIRE 500 ft m Catalog Conte ms 


Column 

Name 

Description 

Null 

1 

Source name 

The identifying natne of the source. A description of source name construction is found above 


2 

R.A. (J2000) 

Right ascension, J2000 (deg) 


3 

Ded. (J200G) 

Declination, J2000 (deg) 


4 

dRA 

Right ascension uncertainty (arcsec ) 


5 

dDecL 

Declination uncertainty (arcsec) 


6 

flux 

Starfinder-extracted flux (mJy) 


7 

une 

Uncertainty of flux (mJy) 


8 

HagFlux 

Flag describing the quality of the flux measurement (see Section 7.3) 


9 

HagStripe 

Flag notes if source is in masked region and may suffer from striping. 1 for masked: 0 for unmasked 


10 

uncSF 

Starlin dei -estimated uncertainty of Flux (mJy) 


n 

sky 

Starfin der-determined sky brightness (MJy sr 1 ) 


12 

fwhm 

Full-widih-at-half-maximum (arcsec) (see Section 7.3) 

0 

13 

flagConf 

Confusion flag (see Section 7). 1 for confused; 0 for not confused 


14 

con- 

Starlin der correlation parameter (see Section 7) 


15 

epochs 

Number of single-epoch images the source was identified in 


16 

s350 

Number of SPIRE 350 ftm Full Lisl sources within 28:4/27” 8 in the LMC'SMC 


17 

ml60 

Number of MIPS 160 ^m Catalog sources wilbin 2 8: 4/27? 8 in Ibc LMC/SMC 



Herschel Space Observatory^ HSO ; II specifies the instrument 
(P for PACS, S for SPIRE) and the wavelength (2, 3 ? 1, 2, 3 
for IO0 S 160, 250, 350, and 500 fim, respectively); PPPP is a 
four-character code for die HERITAGE project, 1IERI; E de- 
notes the epoch (all sources are marked with 0, combined- 
epoch detections), and L designates die list type, 0 for Catalog. 
RR.RRRR and DD.DDD are die right ascension and declina- 
tion, respectively, in degrees to four decimal places, and S is die 
sign of the declination. 
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